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The  purposes  of  this  research  program  v ere  to  extend  or  develop  analytical 
aathoda  for  dot  oral  nine  rotor  blada  aaroalaatic  stability  Halts  and  to 
parfora  stability  celculatlons  ovar  a  range  of  design  and  oparating  vari- 
ablas  for  articulate!  and  nonart lculatad  configurations.  Tha  usafulnass  of 
siaplar  analytical  not  bods  vas  lnvasttgatad  by  coaparlng  results  with  oper¬ 
ating  boundaries  froa  tha  aora  elaborate  analysis. 

In  tbe  part  of  the  investigation  presented  in  this  voluae,  a  staple  energy 
balance  aethod  for  determining  static  stability  vas  used  to  evaluate  rotor 
blade  torsional  divergence  boundar Us,  and  a  set  of  blade  design  charts 
generated.  To  coapare  the  results  of  this  staple  analysis  vith  a  aore  re¬ 
fined  aethod,  and  to  investigate  the  effects  of  parsaeters  not  Included  in 
tbe  staple  aethod,  other  calculations  vere  aade  vith  the  extended  Voraal 
Node  Transient  Analysis. 

One  type  of  practical  rotor  blade  operating  boundary  li  defined  by  peak-to- 
peak  torsional  stress.  The  static  torsional  divergence  stability  boundary 
usually  approximates  this  practical  operating  Halt,  as  verified  by  the 
loraal  Node  Transient  Analysis  calculations.  This  approx iaat ion  is  valid 
for  an  iaportant  nuaber  of  blade  parameters. 

This  part  of  the  investigation  is  Introduced  vith  a  description  of  tbe  tor¬ 
sional  divergence  concept  at  it  applies  to  rotorcraft  blades.  Tbe  condensed 
development  of  the  simple  energy  balance  aethod  is  given.  The  uni  fora  blade 
design  charts  are  tbe  result  of  the  application  of  this  aethod.  The  uniform 
blades  treated  In  the  charts  are  of  arbitrary  sise,  control  systsa  flexibil¬ 
ity,  aerodynamic  root  cutout,  and  section  center-of-gravlty  location.  The 
results  of  additional  calculations  vith  tbe  siaple  aethod  for  nonuniform 
blades  are  also  given  graphically.  The  independent  effects  of  no  nun  1  form 
blade  chord,  torsional  stiffness,  and  inertia  are  shown,  vith  inertial 
effects  reflected  by  the  use  of  tbe  natural  torsional  vibration  aode  as  an 
assumed  divergence  aode.  Tbe  Independent  effects  of  nonuniform  lift-curve 
slope  caused  by  Naeh  number,  pitch-flap  coupling,  and  noosrtl culatlon  are 
presented  In  similar  fashion,  normal  Node  Transient  Analysis  calculations 
for  verification  of  the  simpler  analysis  Include  investigations  of  tbe 
effects  of  nonarticulation,  reduced  air  density,  control  systsa  flexibility, 
aft  center-of-gravlty  position,  and  pitch-flap  coupling  on  articulated 
blades.  The  effects  of  gyro  feedback  coupling  on  nonartlculated  blades  are 
also  presented.  Discussion  of  tbe  results  of  this  part  includes  tbe  signif¬ 
icance  of  the  idealised  torsional  divergence  parsaeters,  the  relative  effects 
on  torsional  divergence  of  the  various  parsaeters,  and  permissible  analyt¬ 
ical  assumptions  useful  in  the  calculation  of  tbe  torsional  divergence 
stability  boundary. 

Recoaaended  preliminary  design  procedures,  incorporating  tbe  charts  gen¬ 
erated  for  this  voluae  are  given  for  the  rapid  definition  of  blades  vith 
adequate  torsional  stlffneso. 
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Th*  investigation  pr««  -nted  la  this  volua*  la  part  of  a  general  study  of 
rotor  blada  aero* lest 1.  Instabilities,  which  la  contained  la  five  rolmm. 
Th*  work  was  performed  under  Contract  DA  AL-177-AJC-3J?(T)  with  th*  U.  8. 
A*my  Aviation  Materiel  Laboratories ,  Fort  Eustls,  Virginia.  Th*  program 
vu  monitored  for  'JSAAVLAAS  by  Mr.  Jo**ph  McOervey. 

Th*  rotor  blad*  torsional  dlv*rg*nc*  analysis,  d*sign  charts,  aad  discus¬ 
sions  in  this  volua*  a-*  th*  result  of  work  don*  at  Sikorsky  Aircraft  by 
Kr.  Charles  f.  Nlebanc*  and  at  th*  Unit*d  Aircraft  Research  Laboratories 
by  Mr.  H.  L.  Claan. 

Volua*  I  of  this  report  contains  th*  development  of  tb*  differential  equa¬ 
tions  of  action  of  an  elastic  rotor  blade  with  cbordwls*  aass  unbalance. 

Volua*  II  presents  a  linearised  discreet  a  sleuth  classical  flutter  analysis 
for  rotor  blade*  with  an  appropriate  parameter  variation  study,  a  comparison 
with  test  data,  aad  a  comparison  with  results  calculated  by  using  the  method 
of  Volua*  I. 

Volt**  III  describes  a  stall  flutter  analysis  based  on  the  calculation  of 
aerodynaalc  work  for  a  cycle  of  blad*  torsional  vibration,  tmo  dimensional 
unsteady  airfoil  test  data  were  used  In  tb*  evaluation  of  tb*  aerodynamic 
work.  Tb*  analysis  was  used  to  jenerat*  stall  flutter  boundaries. 

Volume  V  presents  tb*  results  of  a  study  of  flapping  and  coupled  flap- lag 
1  ratability.  The  results  of  a  parametric  study  based  on  a  single  degree  of- 
freedra  flapping  or  flatwise  bendiug  analysis  are  presented.  Comparisons 
are  aad*  with  results  from  th*  aore  elaborate  method  of  Volua*  I.  Tb*  re¬ 
sults  obtained  by  using  the  method  of  Volume  I  to  determine  the  coupled 
flap-lag  response  of  a  rotor  to  a  number  of  sudden  control  changes  are  pre¬ 
sented. 


v 


TABLE  OF  OOWEWD 


SliMARY .  Ill 

FOREWORD  .  v 

LIST  OP  ILLUSTRATIONS . rill 

LIST  OF  TABLES .  *1 

LIST  OP  SYMBOLS . xli 

INTRODUCTION .  1 

DEVELOPMENT  AND  DISCUSSION  OP  TORSIONAL 

DIVERGENCE  ANALYSIS .  3 

TORSIONAL  DIVERGENCE  DESIGN  CHARTS .  23 

TORSIONAL  DIVERGENCE  CALCULATIONS  WITH  THE  EXTENDED 

NORMAL  NODE  TRANSIENT  ANALYSIS .  31  . 

DISCUSSIONS  OP  TORSIONAL  DIVERGENCE  RESULTS .  38 

CONCLUSIONS . 108 

RECOMMENDATION . 109 

DISTRIBUTION .  110 


v  I  i 


LI8T  or  ILLUSTRATIONS 


TlW  PSfie 

1  Comparison  of  Differential  Equation  Solution 

and  Assumed  Divergence  Mode . 1*7 

2  Advance  Ratio  Versus  Torsional  Stiffness  Coefficient 

for  Various  Control  Stiffness  Ratios  and  Center-of- 
Oravity  Positions;  X0s  2  0  ,  o0  =  2 rr . 1*8 

3  Reciprocal  of  Control  Stiffness  Ratio  Versus 

Torsional  Stiffness  Coefficient;  X0$s  0  ,  o0  *  2 ir  ....  6U 

U  Center-of-Gravity  Position  With  Respect  to  75% 

Chord  Versus  Torsional  Stiffness  Coefficient; 

X©$  *  0  ,  o0  s  21T . 69 

3  Effect  of  Linear  Variation  in  Chord  on  Retreating 

Blade  Torsional  Divergence;  l/K  =  .001  ,  Yc0/C  =  0  , 

X0*  *  .25  ,  Xog  «  0  ,  o0  s  2tr . 80 

6  Effect  of  Structural  Cutout  on  Retreating  Blade 
Torsional  Divergence;  l/K  *  .001  ,  Ycq/C  =  0  , 

Xoa  *  0  ,  o0  *  2ir . 8l 

7  Effect  of  Linear  Variation  in  Stiffness  on 
Retreating  Blade  Torsional  Divergence; 

l/K  *  .001  ,  Yco/C  *  0  ,  X0*  =  0  ,  X0$  =0  ,  o0  «  Ztt . 82 

8  Effect  of  Inertial  Cutout  on  Retreating  Blade 
Torsional  Divergence;  l/K  *  .001  ,  Ycq/C  *  0  , 

X0s  s  0  .  °o  s  2ir . 83 

9  Effect  of  Linear  Variation  in  Inertia  on 
Retreating  Blade  Torsional  Divert* e ice  ; 

l/K  =  .001 ,  Xo*  =  0  ,  XM  =  0  ,  a0  =  27r . 81* 

10  Effect  of  Prandtl-Glauert  Mach  Number  Correction 
on  Advancing  Blade  Torsional  Divergence; 

l/K  5  .001  ,  Ycq/C  s  -  .10  ,  Xqa  =  0  ,  Xqj  =  0  •  •  •  . 85 

11  Effect  of  Prandtl-Glauert  Mach  Number  Correction 
on  Retreating  Blade  Torsional  Divergence; 

l/K  =  .001  ,  Yco/C  =  0  ,  X0A  =  0  ,  XP#  =  0  . 86 

12  Retreating  Blade  Torsional  Divergence  on 
Articulated  and  Nonart iculated  Blades; l/K =.001, 

Yco/C  =  .005  ,  Xw  =  0  ,  Xos  =  0  ,  ^OoCoYm  =  .  168 . 87 

13  Advancing  Blade  Torsional  Divergence  on  Articulated 
and  Nonarticulated  Blades;  l/K  =.001  ,  Ycq/C  =  -.005, 

X  qi  *  0  ,  Xqj  -  0  i  p  o0C0  /m  ~  .  1 68  ••••••••.•••  88 

viii 


Figure  Pm 


Ik  Effect  of  Pitch-Flap  Coupling  on  Retreating 

Blade  Torsional  Divergence;  l/K  =  .001  , 

X<>a  -  0  ,  Xqs  5  0  ,  o0  5  2tt . 89 

15  Blade  First  Torsion  Modal  Response  Versus  Azimuth 
Angle  for  the  Basic  Articulated  Rotor;  pz  1.6, 

0.76*  =0“  ,  Ycg/C  =  0  ,  M|(t0  =  85 . 90 

16  Blade  First  Torsion  Modal  Response  Versus  Azimuth 
Angle  for  the  Articulated  Rotor  with  Aft  Center 
of  Gravity;  p  -  .6  ,  07#*  =  2°  ,  Yc6/C  =  -.15  , 

M,,to  =  .85 . 91 

17  Vibratory  Stress  Versus  Advance  Ratio  for  the 
Basic  Articulated  Rotor;  l/K  =  0,  Ygq/C  =  0, 

X0a  =  .12  ,  X0,  =  .060  ,  M,ti0  =  .85  ,  8,  =  0  ,  Ooft'S*  =  96 . 92 

18  Vibratory  Stress  Versus  Advance  Ratio  for  the 
Basic  Articulated  Rotor  with  Reduced  Air 
Density;  l/K  =  0  ,  YGG/C  =  0  ,  X0A  =  .  12  , 

Xos  =  060  ,  M,  *  =  .85  ,  p  =  .001585  SLUGS/ FT3  , 

83  =  0°  ,  a0n*S*  =  144 . 93 

19  Vibratory  Stress  Versus  Advance  Ratio  for  the 

Basic  Articulated  Rotor  With  Control  System 
Flexibility;  l/K  =  .2  ,  YCg/C  =  0  ,  =  .12  , 

X0#  =  .060  ,  M|,»0  =  .85  ,  8,  =0°  ,  OoQ'S,,  =  96 . 9k 

20  Vibratory  Stress  Versus  Advance  Ratio  for  the 
Basic  Articulated  Rotor  With  Pitch-Flap  Coup¬ 
ling;  l/K  =0  ,  Ycg/C  =0  ,  X0*  =  .12  ,  X0,  =  .060  , 

Ml|t0  =  -85  ,  8,  =  45°  ,  a0a*S*  =  96 . 95 

21  Vibratory  Stress  Versus  Center -of -Gravity 
Offset  for  the  Articulated  Rotor;  l/K  =  0  , 

Xo»=  .12  ,  Xm  =  .060,  p.  -  .6  ,  Mlt#0  =  -85  , 

S3  =  0  ,  aoa*S*  =  96  ,  0.75R=  2° . 96 

22  Vibratory  Stress  Versus  Advance  Ratio  for  the 
Basic  Nonaxticulated  Rotor;  l/K  =0  ,  YGG/C  =  0  , 

X0a  =  .12  ,  Xcs  =  .060  ,  M,|V0  =  .85  ,  oea*S*  =  96  . 97 

23  Vibratory  Stress  Versus  Center-of-Gravity  Offset 
for  the  Nonarticulated  Rotor;  l/K  =  0  ,  X0A  =  .12  , 

Xos  =  .060  ,  p  =  .6  ,  MItf0  =  .85  ,  ooa2S*  =  96  , 

0.76*  =  2°  ,  A,s  =  -1°  ,  Bis  =  2 . 98 


ix 


Figure 

2k 


Blade  First  Torsion  Modal  Response  Versus  Azimuth 
Angle  for  the  Nonarticulated  Rotor;  Ycq/C  =  0  , 
fj,  s  \  . A  ,  M|(f0  s  . 85 . 


Page 


i 


99 


25  Effects  of  Torsional  Stiffness  Coefficient 
Variations  on  Typical  Torsional  Divergence 
Flight  Condition  Boundaries;  l/K  =  .001  ,  Ycq/C  =  0  , 

X0A  =  .25  ,  X0,  =  0  ,  o0  =  27 r  ,  S„(Ra)*  =  2GJ/ pa0CfR2  •  •  •  •  100 

26  Effects  of  Chordvise  Center -of -Gravity 
Location  on  Typical  Torsional  Divergence 
Flight  Condition  Boundaries;  l/K  ~  .001  .  X0a  =  .25  , 

Xos  s  0  ,  c0  =  2t r  ,  S^Rfl)1  =  2GJ/>o0Co*R4  =  16 .750FT/SEC*  .  .  101 


27  Effects  of  Aerodynamic  Root  Cutout  on  Typical 
Torsional  Divergence  Flight  Condition 
Boundaries;  l/K  =  .001  ,  Ycq/C  =  0  ,  XM  =  0  ,  o0  =  2tt  , 

Sft(Ra)*  *  2Qyi/po0C^Rl  -  16,750  FT /SEC* . 102 

28  Effects  of  Control  Stiffness  on  Typical 
Torsional  Divergence  Flight  Condition 
Boundaries;  Yco/C  =  0  ,  X0A  =  .25  ,  XM  =  0  ,  o0  -  2tt  , 

Sn(Rfl)*  *  2GJ//>  OoCo*  R*  *  16,750  FT/ SEC* . 103 


> 

I 


i 

■j 


x 


{ 

1 


LIST  OF  TABLES 


Table  Paste 


I  Basic  Properties  of  Rotor  System  Used  in  Torsional 

Divergence  Studies  With  the  Normal  Mode  Transient 
Analysis .  lOfc 

II  Rotor  Configurations  for  Torsional  Divergence 

Study  With  the  Normal  Mode  Transient  Analysis .  105 

III  Flight  Conditions  for  Torsional  Divergence 

Study  With  Normal  Mode  Transient  Analysis .  106 

IV  Configuration  and  Flight  Condition  Combinations 

for  Torsional  Divergence  Study  With  the  Normal 
Mode  Transient  Analysis  .  .  .  . .  107 


xi 


ipttd  of  sound,  ft/sec 
lateral  cyclic  pitch,  deg 


A* 

Ai» 

Oo  reference  value  (at  75S  radius)  of  slope  of 

C  versus  Qf  curve,  dC^/ddr  rad-1 

o  local  value  of  slope  of  Cj  versus  dr  curve, 

dCj/ddr  rad-1 

Bis  longitudinal  cyclic  pitch,  deg 

Cues  blade  nodal  constants  defined  in  Eq.  (W  through  Eq.  (53) 

C  local  value  of  blade  chord,  ft. 

Co  reference  value  (at  75*  radius)  of  blade  chord,  ft. 

C  local  chord  divided  by  rotor  radius 

Cg  section  drag  coefficient,  d/{'/ipU*C) 

Cmc/i  section  pitching  nonent  coefficient  about  the 

25*  chord ,  me/4  /  ( Wz p UzCf ) 

C f,  section  lift  coefficient,  //(|/2/»U2C) 

E  Young's  Modulus  of  elasticity,  lb/ft2 

Elf  flatwise  blade  bending  stiffness,  lb-ft2 

So  reference  distance  (at  75*  radius)  between  the  coin¬ 

cident  elastic  and  center-of-gravity  axes  and  the 
aerodynamic  center,  ft 

V  flapping  binge  radius  divided  by  rotor  radius 

e«  distance  of  the  centroid  of  blade  spar  area  forward 

of  the  elastic  axis,  ft. 

p 

G  shear  modulus  of  elasticity, lb/ft 

p 

CJ  local  blade  torsional  stiffness,  lb-ft 

GJt  local  blade  torsional  stiffness  at  tip,  lb-ft2 

GJft  local  blade  torsional  stiffness  at  root,  lb-ft2 


xil 


G  Jo 


ff3o 


g 

r. 

J 

K 

K, 

k 

*r 


Le 

i 

Mf 

M 

Mi.*o 
M  »*.*o 
M  n.tro 
Me 

m 


iflj 


fh 

m4 


reference  (at  75 S  radiua)  blade  torsional 
stiffness,  lb-ft2 

nond'aensional  blade  torsional  stiffness , 
defined  in  Eq.  (98) 

acceleration  of  gravity,  ft/sec2 

nondioensional  blade  flapping  moment  of  inertia 
defined  in  Eq.  (99) 

I, 

torsional  stiffness  constant  of  blade  section,  ft 

control  stif fness-to-blade  stiffness  ratio, 

K  --  K,  lb  /GJ0 

control  stiffness  in  torsion,  ft-lb/rad 

radius  of  gyration  of  blade  section  aass,  ft. 

radius  of  gyration  of  spar  cross-sectional  area,  ft. 

local  blade  Bass  radlur.  of  gyration  in  torsion,  ft. 

length  of  blade  outboard  of  push  rod  location,  ft. 

lift  force  increment  divided  by  blade  radial 
lncreeent,  lb/ft 

blade  eass  divided  by  m<,R 

blade  Bass  ratio  Mr  s  "'o/pir  C 

nondlmens ional  cooent 

advancing  blade  tip  Mach  nuaber 

advancing  blade  Mach  nuaber  at  the  75 t  rotor  radius 

retreating  blade  Mach  nuaber  at  the  75S  rotor  radius 

noaent  due  to  elastic  deformation  of  the  blade,  ft-lb 

blade  Bass  increment  divided  by  radial  increment, 
slugs/ft 

reference  value  of  m  ,  olugs/ft 
local  value  of  m  divided  by  m0 

section  aerodynanlc  pitching  moment  per  unit  span  about 
elastic  axis,  ft-lb/ft 


xiil 


mc/4  steady  state  section  pitching  moment  per  unit  span,  ft-lb/ft 

"km  counterweight  mass  per  unit  span,  lb-sec2/ft2 

Po  •  Pi  dlaensloniess  constants  in  Eq.  ( IO9)  which  are  determined 
by  the  boundary  conditions  for  the  differential  equation 

Q  local  dynaalc  pressure,  lb/ft2 

Of,  first  torsion  modal  amplitude,  rad 

qV|  first  flatwise  bending  aodal  amplitude,  ft 

qy(  )  amplitude  of  (  )th  edgewise  deflection  mode  (equals  blade 

tip  deflection  in  direction  divided  by  R  when 

1  at  tip  is  defined  as  1.0} 

qw(  )  amplitude  of  (  )th  flatwise  deflect' ->n  mode  (equals  blade 
tip  deflection  in  1,  direction  divided  by  R  when 
y«v(  >  at  tip  is  defined  as  1.0) 

qf(  )  amplitude  of  (  )th  elastic  torsional  mode  (equals  elastic 
twist  angle  about  >•  axis  in  radians  at  tip  when  Yg{  , 
at  tip  is  defined  as  1.0) 

R  rotor  radius,  ft 

Rj  tip  value  of  blade  torrional  Inertia  per  unit  span  divided 

by  root  value 

r  distance  along  the  blade  radius  from  the  flapping  hinge 

for  the  articulated  blades  and  from  the  rotor  axis  for 
the  nonart lculated  blade,  ft 

F  nondlmensional  value  of  r  i  P  *  r/R 

rT  value  of  r  at  blade  tip,  ft 

r0Cw  value  of  t  at  inboard  end  of  counterweight,  ft 

FCg  radius  of  blade  center  of  gravity  defined  in  Eq.  (100) 

S*  aerodynamic  shear  force  per  unit  span,  lb/ ft 

St  nondlmensional  parameter,  Sc  *  GJo/p^c! R2( A*)2 

SR  nondlmensional  torsional  stiffness  coefficient, 

defined  in  Eq.  (106) 

t  time,  sec 

U  resultant  velocity  at  blade  section,  U*  •/ \jl  +  U7 

xl  V 


, ft/sec 


Up  velocity  component  at  blade  section,  normal  to  the  local 

elastic  axis  direction  and  Ut  *  ft/sec 

UT  velocity  component  at  blade  section,  noraal  to  tbe  local 

elastic  axis  direction  and  parallel  to  tbe  plane  of  rota¬ 
tion,  ft /sec 

Of  local  nondlaenslonal  relative  air  velocity,  noraal  to  tbe 

undeflected  blade  radius  and  parallel  to  tbe  plane  of 
rotation.  Of*  I  ♦  r  ♦fism* 

V  aircraft  forward  velocity,  knots  or  ft/sec 

V#  local  edgewise  elastic  dlsplaceaent  of  elastic  axis,  ft 

We  local  flatwise  elastic  dlsplaceaent  of  elastic  axis,  ft 

■(uVa,**  coordinates  In  the  nth  axis  systea 

i  radial  distance  froa  rotor  axis  of  rotation  divided  by 

rotor  radius 

Xqa  radius  of  innermost  aerodynaalcally  effective  blade 

section,  divided  by  rotor  radius  (aerodynamic  root  cutout) 

Xom  radius  of  lnneraost  blade  section  with  appreciable  torsional 

inertia,  divided  by  rotor  radius  (inertial  root  cutout) 

Xqs  radius  of  lnneraost  torslonally  flexible  blade  section, 

divided  by  rotor  radius  (structural  root  cutout) 

*ac  local  chordvlse  distance  frca  the  elastic  axis  to  tbe 

aerodynamic  center  divided  by  rotor  radius,  positive 
when  the  aerodynamic  center  Is  forward 

Yco  local  chordvlse  distance  from  blade  elastic  axis  to  blade 

center  of  gravity  position,  positive  when  center  of  gravity 
Is  forward,  ft 

Ygc  nondlaenslonal  center  of  gravity  offset,  Yco/R 

a,  rotor  shaft  angle  of  attack,  positive  nose  up,  rad 

Of  local  blade  section  angle  of  attack,  rad 

0  blade  flapping  angle  at  hinge,  rad 

fy,  first  uncoupled  torsional  natural  vibration  normalized 

aodc  shape 

XW|  first  uncoupled  bending  natural  vibration  normalized 

node  shape 

Y$ ,  first  derivative  of  Yq,  ,  ^Yg,  /& 
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first  derivative  of  Ym,  ,  dyW|/dV 

value  of  fii,  at  push  rod  location 

value  of  at  push  rod  location 

(  )th  uncoupled  edgewise  natural  vibration 
norma  Used  mode  shape 

(  )th  uncoupled  flatwise  natural  vibration 
normalised  mole  shape 

(  )th  uncoupled  torsional  natural  vibration 
normalised  node  shape 

distance  between  chordwise  location  of  counterweight  and 
chordvlse  location  of  spar  centroid,  positive  when  former 
is  between  centroid  and  leading  edge,  ft 

blade  lead  angle,  rad 

pitch-flap  coupling  angle,  S»  *  orcton  dQ/d$ 
total  local  blade  pitch  angle,  rad 

local  blade  pitch  angle  due  to  control  system  Input  and 
built-in  twist,  rad 

local  blade  pitch  angle  due  to  built-in  linear  twist,  rad 
local  elastic  twist  angle,  rad 

rate  of  change  of  0%  with  respect  to  1  ,  rad 

collective  pitch  setting  expressed  at  the  75f  radius,  deg 

total  air  velocity  component  parallel  to  rotor  shaft 
divided  by  rotational  tip  speed  HR 

advance  ratio,  V/flR 

dummy  variable  of  integration  along  rigid  blade  axis,  ft 
air  density,  alugs/ft^ 
local  Inflow  angle,  rad 

blade  aslmuth  angle,  positive  in  direction  of  rotation 
from  downwind,  rod  or  deg 

first  torsional  natural  frequency  divided  by  rotor 
rotational  frequency 


Qtl  first  flatwise  banding  natural  frequency  divided  by  rotor 

rotational  frequency 

u>v,  ,  uncoupled  natural  frequency  of  (  )lh  edgewise  vibratory 

■ode,  rtd/sec 

u/,(  ,  uncoupled  neturel  frequency  of  (  )th  flatwise  vlbretory 

■ode,  rod/ sec 

utg{  |  uncoupled  neturel  frequency  of  (  )*h  torelonel  vlbretory 

■ode,  red/aec 

0  rotor  angular  velocity,  red/aec 

SUBSCRIPTS 

i**0  indice tea  axis  system  under  cons ideret Ion,  es  defined 

In  Volume  I  of  this  report 

C/a  indicates  quantity  evaluated  at  2$%  chord 

K/e  indicates  quantity  evaluated  at  75$  chord 

a  denotes  a  known  discreet  set  of  values,  as  used  In  Eq.  (110) 

eg  Indicates  quantity  evaluated  at  local  chordvlse  center  of 

gravity 

c«  indicates  quantity  pertaining  to  the  blade  counterweight 

eo  indicates  quantity  evaluated  at  elastic  axis 

ro  indicates  flap  damper  or  quantity  evaluated  at  flap 

damper  attachment  point  as  appropriate 

•  flatwise  node  subscripts 

I  .I'.h  torsional  node  subscripts 

lO  indicates  lag  damper  or  quantity  evaluated  at  lag 

damper  attachment  point  as  appropriate 

p*  indicates  push rod  or  quantity  evaluated  at  pushrod 

attachment  point  as  appropriate 

p,p',n  edgewise  mode  subscripts  (maximum  value  ■  2) 

t  !  0  indicates  quantity  evaluated  at  t  •  0 


■,V ,*  Indicates  quantity  evaluated  in  »  ,  y  or  t 

direction  as  defined  in  Volume  1  of  this  report 
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*()•*(**()  indicates  quantity  evaluated  In  »  ,  y  or  i 

direction  of  (  )  axis  syntoi  a*  defined  In  Voluae  I  of 
this  report 

B  ®  indicates  quantity  evaluated  at  yM(f  *  0 
DOIVATIVI  MOT  ATI  OB 


d(  )/dr 

indicataa  derivative  of  ( 

)  vith  respect  to  r 

i  >' 

indicataa  derivative  of  ( 

)  with  respect  to  F 

(') 

indicates  derivative  of  ( 

)  with  respect  to  t 

(*) 

indicates  derivative  of  ( 

)  with  respect  to  i 

(l.s.,  t  ) 

inOTT.umi|8 


n 

indicates 

n 

Indicates 

of  factors 

vac  tor  quantity 

quantity  nondiXMnsiooalised  through  tba  usa 
,  R  ,  m0  ,  and  A 
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INTRODUCTION 


One  of  the  basic  aeroclastic  investigations  applying  to  fixed  wings  is  the 
torsional  divergence  analysis,  for  cases  where  the  center  of  pressure  is 
ahead  of  the  elastic  axis.  The  central  concept  is  the  presence  of  an  aero¬ 
dynes  ic  torsional  couple  on  the  wing,  which  increases  linearly  with  suoglm 
of  attack  and  which  therefore  varies  in  similar  fashion  as  the  wing  twints. 
These  couples,  which  are  opposed  by  an  elastic  restoring  torque.  Increase 
with  the  square  of  the  forward  velocity  for  a  given  angle  of  attack.  In 
such  a  case,  the  aerodynaaic  couple  is  proportional  to  the  square  of  the 
speed  and  linear  with  the  twisting  deformation,  while  elastic  resisting 
couple  is  proportional  to  the  twisting  deformation  alone.  It  is  clear 
that  a  velocity  will  exist  for  which  an  Increment  in  the  aerodynaaic  twist¬ 
ing  couple  is  equal  to  the  corresponding  increment  of  elastic  resisting 
couple  for  any  arbitrary  deflection.  This  is  known  as  the  torsional  diver¬ 
gence  speed.  Below  this  speed,  the  elastic  resistance  to  an  increment  in 
torsional  deflection  will  be  greater  than  the  corresponding  aerodynamic 
couple,  and  the  deflection  will  find  its  equilibrium  value.  Above  it,  the 
Increment  in  elastic  resistance  will  be  less  than  the  increment  in  the 
aerodynamic  couple,  and  the  simple  divergence  theory  predicts  a  continuous 
growth  in  deflection.  It  should  be  noted  in  passing  that  the  divergence 
speed  may  be  either  larger  or  smaller  than  the  classical  flutter  speed. 

The  above  concept  of  a  torsional  divergence  for  a  fixed  wing  led  to  an 
examination  of  similar  situations  existing  for  helicopter  rotors.  Obvioutfy, 
the  velocity  distribution  on  a  helicopter  rotor  varies  along  the  span  and 
is  rapidly  and  continuously  changing.  Hence,  the  static  stability  analysis 
for  torsional  divergence  applies  only  to  an  Instantaneous  condition.  In 
fact,  the  torsional  divergence  situation  for  a  helicopter  rotor  blade 
usually  develops  on  the  retreating  blade  for  advance  ratios  greater  than 
unity.  The  blade  is  then  traveling  backwards  (sharp  edge  first)  through 
the  air  for  part  of  each  revolution,  and  the  aerodynamic  center  moves  close 
to  what  is  normally  the  75%  chord  position.  This  produces  a  large  tor¬ 
sional  moment  arm  about  the  blade  elastic  axis  and  center-of-gravity  posi¬ 
tion  at  or  near  the  25%  chord  position.  Hence,  torsional  divergence  will 
be  encountered  for  the  retreating  blade  even  though  the  relative  velocity 
is  comparatively  low. 

A  more  detailed  description  of  the  theoretical  approach  to  the  torsional 
divergence  calculations  will  follow  later.  However,  a  definition  of  tor¬ 
sional  divergence  can  be  given  at  this  point.  In  terms  of  theory,  the 
torsional  divergence  boundary  is  that  set  of  flight  conditions  where  blade 
static  aeroelastlc  stability  in  torsion  is  Just  satisfied  at  either  the 
advancing  or  retreating  blade  positions.  The  practical  effect  of  exceeding 
this  boundary  has  been  found  to  be  a  steep  rise  in  blade  torsional  vibra¬ 
tory  stress  levels. 


1 


In  this  volume  of  the  report,  the  flight  conditions  for  torsional  diver¬ 
gence  have  been  examined  extensively,  and  useful  design  charts  have  been 
prepared.  These  critical  flight  conditions  have  also  been  examined  vith 
the  much  more  rigorous  extended  Normal  Mode  Transient  Analysis,  vhich 
provides  solutions  for  the  flexible  blade  equations  of  motion.  It  is  gen¬ 
erally  verified  that  the  torsional  divergence  analysis  fulfills  its  pur¬ 
pose,  vhich  is  the  definition  of  potentially  troublesome  flight  conditions 
which  should  be  avoided  or  investigated  further. 
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DEVELOPMENT  AND  DISCUSSION  OF  TORSIONAL  DIVERGENCE  ANALYSIS 


Rotor  blades  on  high  speed  helicopters  can  experience  torsional  divergence 
as  discussed  earlier.  One  essential  difference  between  helicopter  blade 
and  fixed  wing  torsional  divergence  is  the  role  played  by  the  center-of- 
gravity  location.  For  a  fixed  wing,  the  torsional  moment  causing  the 
divergence  is  the  lift  force  increment  multiplied  by  the  distance  between 
the  aerodynamic  center  of  pressure  and  the  elastic  axis.  For  articulated 
rotor  blades,  and  nonarticulated  rotor  blades  of  conventional  stiffness, 
most  of  the  resistance  to  blade  vertical  deflection  is  supplied  by  centri¬ 
fugal  force  component?  rather  than  by  elastic  forces.  For  this  reason, 
the  torsional  moment  causing  the  torsional  divergence  of  a  rotor  blade  of 
conventional  stiffness  is  very  nearly  the  lift  force  increment  multiplied 
by  the  distance  between  the  aerodynamic  center  of  pressure  and  the  center 
of  gravity.  On  a  blade  of  conventional  stiffness,  the  distance  between 
the  elastic  axis  and  the  aerodynamic  center  will  have  very  little  effect 
on  the  torsional  divergence  speed. 

In  the  case  of  nonarticulated  rotor  blades  which  are  to  be  stopped  in 
flight,  the  torsional  divergence  mechanism  is  more  like  that  for  fixed 
wings.  As  the  rotor  angular  velocity,  is  decreased,  the  blade  vertical 
stiffness  due  to  centrifugal  effects  diminishes,  and  the  vertical  elastic 
forces  must  resist  a  larger  portion  of  the  vertical  loads  on  the  blade. 

If  the  rotor  is  stopped,  the  basic  torsional  divergence  mechanism  is  the 
same  as  that  for  a  fixed  wing.  It  should  be  noted  that  the  combined 
effects  of  centrifugal  and  elastic  stiffness  are  considered  in  the  non¬ 
articulated  blade  torsional  divergence  analysis  to  be  presented  later. 

As  explained  above,  the  torsional  moment  which  causes  blade  torsional 
divergence  arises  because  the  lift  forces  on  the  blade  are  not  coincident 
with  the  center  of  centrifugal  and  elastic  forces  resisting  the  vertical 
deflection  of  the  blade.  In  order  for  the  moment  to  produce  a  potentially 
unstable  mechanism,  the  moment  must  cause  the  blade  angle  of  attack  to 
become  greater  when  the  lift  force  becomes  greater.  This  condition  is 
obtained  when  the  aerodynamic  center  of  pressure  is  forward  cf  the  center 
of  the  effective  centrifugal  and  elastic  forces  resisting  vertical  motion 
of  the  blade.  If  the  aerodynamic  center  of  pressure  is  aft  of  the  center 
of  effective  vertical  resistance,  torsional  divergence  cannot  occur,  since 
an  increase  in  lift  will  result  in  a  blade  torsional  deflection  which  will 
decrease  the  angle  of  attack. 

It  is  instructive  to  examine  the  conditions  for  which  a  conventional  blade 
can  encounter  the  stable  and  unstable  situations  Just  described.  If  the 
blade  center  of  gravity,  aerodynamic  center  of  pressure,  and  elastic  axis 
are  at  the  25%  chord  position,  ro  torsional  divergence  is  possible  on  an 
advancing  blade,  since  the  aerodynamic  center  of  pressure  coincides  with 
the  center  of  vertical  resistance.  If  the  center  of  gravity  is  ahead  of 
the  25%  chord  position,  the  moment  produced  by  a  lift  increment  is  stat¬ 
ically  stabilizing.  In  this  case,  the  blade  on  the  advancing  side  of  the 
rotor  will  tend  to  twist  so  as  to  reduce  the  lift  increment  due  to  a  gust 
or  a  control  movement.  If,  on  the  other  hand,  the  center  of  gravity  is 
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aft  of  the  25%  chord  position,  the  moment  produced  by  a  lift  increment  is 
potentially  unstable  for  the  advancing  blade.  As  pointed  out  earlier,  the 
relative  position  of  the  elastic  axis  with  respect  to  the  aerodynamic 
center  has  a  similar  but  much  weaker  effect  for  the  conventional  rotor 
blade.  Thus  the  usual  rotor  blade  configuration  will  be  Immune  to  tor¬ 
sional  divergence  on  the  advancing  side  of  the  rotor.  This  is  not  true, 
however,  for  a  blade  on  the  retreating  side  of  the  rotor  when  the  advance 
ratio  is  high  enough  to  produce  reverse  flow  over  most  of  the  blade.  In 
the  reverse  flov  region,  the  blade  is  traveling  through  the  air  backwards, 
and  the  aerodynamic  center  of  pressure  moves  to  the  position  vhich  is 
usually  referred  to  as  the  75%  chord  position.  Since  the  blade  center  of 
gravity  and  elastic  axis  remain  at  the  25%  chord  position,  the  distance 
between  the  aerodynamic  center  of  pressure  and  the  center  of  the  effective 
blade  vertical  stiffness  is  50%  of  the  chord,  This  produces  a  potential 
static  Instability,  which  fortunately  is  mitigated  by  the  lower  dynamic 
pressures  on  the  retreating  blade.  It  should  be  noted  that  for  the  re¬ 
treating  blade,  moving  the  center  of  gravity  towards  the  normal  trailing 
edge  decreases  the  unstable  tendency,  since  it  is  then  closer  to  the  actual 
aerodynamic  center.  The  opposite  effect  is  produced  by  moving  the  center 
of  gravity  forward  towards  the  normal  leading  edge. 

As  mentioned  in  the  previous  discussion,  the  usual  rotor  blade  configura¬ 
tion  will  experience  static  torsional  divergence  only  in  a  region  of  re¬ 
verse  flow.  A  single  degree  of  freedom  analysis  is  adequate  for  a  study 
of  the  effect  of  torsional  stiffness  on  torsional  divegence  of  an  artic¬ 
ulated  blade  with  this  normal  configuration.  In  order  to  evaluate  the 
effects  of  pitch-flap  coupling  and  nonarticulation  with  arbitrary  center 
of  gravity  location,  a  two-degree-of- freedom  analysis  is  required,  and 
was  therefore  used  in  this  investigation. 

The  general  approach  followed  in  this  analysis  is  the  development  of  two 
equations  of  static  stability  involving  one  torsional  and  one  flapping, 
or  flatwise  bending,  degree  of  freedom.  The  determinant  of  these  two 
homogeneous  equations  is  evaluated.  A  zero  value  for  this  determinant 
indicates  that  the  blade  is  operating  on  the  boundary  between  a  statically 
stable  and  a  statically  unstable  regime.  Operatior  of  the  rotor  in  a 
flight  condition  beyond  this  static  stability  boundary  does  not  imply  that 
blade  deflections  will  increase  indefinitely,  since  the  condition  will 
generally  exist  over  only  a  portion  of  the  rotor  disc.  The  amplitude  of 
the  blade  deflections  is  affected  not  only  by  the  extent  and  severity  of 
the  statically  unstable  region,  but  also  by  blade  inertia,  damping,  and 
blade  motions  upon  entry  into  the  unstable  region.  Therefore,  the  results 
of  an  analysis  of  the  type  presented  here  serve  principally  to  indicate 
those  flight  conditions  which  are  potentially  troublesome  from  the  stand¬ 
point  of  adequate  blade  torsional  stiffness  end  which  should  be  analyzed 
in  further  detail. 

Inclusion  of  chordwise  mass  unbalance  for  nonarticulated  blades  and  pitch- 
flap  coupling  effects  requires  consideration  of  blade  deflection  out  of 
the  plane  of  rotation  as  well  as  blade  torsional  deflections  about  the 
elastic  axis.  It  was  assumed  that  the  first  natural  torsional  vibration 
mode  and  the  rigid  flapping  mode  would  adequately  represent  the  deflections 
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of  articulated  rotors.  It  vat  further  assumed  that  the  first  torsion  and 
flatwise  bending  modes  would  repreeent  the  deflection  of  nonart lculat ad 
rotors.  The  total  torsional  spring  coefficient  was  obtained  from  the 
blade  equations  of  motion  presented  in  Volume  I  of  this  report.  A  list  of 
the  principal  simplifying  assumptions  made  to  arrive  at  the  equations  used 
here  is  given  below: 

1.  All  dynamic  terms  are  neglected. 

2.  All  edgewise  and  lag  terms  are  neglected,  as  are  the  terns 
Involving  flatwise  and  torsional  modal  amplitudes  other 
than  the  first. 

3.  For  articulated  rotors,  the  flatwise  bending  amplitude 
is  zero;  for  nonart iculated  rotors,  the  blade  is  rigidly 
cantilevered  at  the  root. 

The  lift  force  is  the  only  aerodynamic  force.  It  is 
described  by  means  of  a  linearized  lift-curve  slope, 
which  is  the  same  absolute  magnitude  in  both  conven¬ 
tional  and  reverse  incompressible  flows.  When  desired, 
compressibility  effects  can  be  accounted  for  by  the 
Prandtl-Glauert  correction,  which  is  applied  to  both 
conventional  and  reverse  flow. 

5-  Only  steady-state  aerodynamics  are  considered  and  the 
aerodynamic  center  of  pressure  is  located  at  the  25% 
chord  for  conventional  flow  or  at  the  15%  chord  for 
reverse  flow. 

6.  The  resultant  section  velocity  is  equal  to  the  tangential 
velocity  component. 

7.  Terms  involving  collective  or  cyclic  pitch,  blade  pre-twist, 
preconing  and  second  and  higher  products  of  modal  amplitudes 
are  neglected. 

8.  Blade  angles  of  attack  remain  small  in  forward  flow  and  are 
close  to  180°  in  reverse  flow. 

The  equations  of  static  stability  were  developed  by  applying  the  assump¬ 
tions  mentioned  above  to  modal  equations,  used  in  the  Normal  Mode  Transient 
Analysis  of  Volume  I.  These  modal  equations  arise  by  effectively  setting 
the  virtual  work  done  by  the  aerodynamic  loads  and  centrifugal  effects 
equal  to  the  change  in  strain  energy  during  a  virtual  displacement  of  the 
assumed  modes.  This  consideration  of  energy  terms  gives  rise  to  the 
description  of  the  static  stability  analysis  as  an  "energy  method". 


The  development  of  the  analysis  will  begin  with  the  modal  equation*  givej 
In  Voluae  X.  The  equation*  needed  are  the  flatwise  bending  nodal  equation, 
the  torsional  modal  equation,  and  the  flap  angle  nodal  equation.  In  Volume 
X  these  appear  as  Eqs.  (80),  (89),  and  (91)  respectively.  For  convenience, 
they  are  reproduced  below  as  Eqs.  (1),  (2),  and  (3).  Many  of  the  terms  In 
these  equations  from  Volume  I  ore  eliminated  herein,  since  they  do  not  per¬ 
tain  to  the  static  divergence  concept.  Some  terms  will  be  used  fur  the 
divergence  analysis  of  both  articulated  and  nonartlculated  blades,  while 
others  will  be  applicable  to  articulated  blades  only,  or  to  nonartlculated 
blades  only.  The  specific  reasons  for  deleting  the  various  terms  will  be 
given.  Volume  I,  of  course,  gives  a  more  complete  discussion  of  the  modal 
equations. 

The  flatwise  bending  modal  equation  from  Volume  I  of  this  report  is 


0  ’  /*.,  [iS»>,,co*e-l5,l>5«ine]  df-Mp,co«8M(r.')w*«'n9Lo“tol)'4.. 

0 

♦  tint  ■  eo.M,)>] 

-  ♦  z  ll-CO«2$,)--±-(l+2»)«in2q] 

-«,co.2«(,(c5-or5CJ)|<-c„^+  28,(1  ^*)ci|t-e,c,it} 

♦Z  iV([-2t,ne0C,4^-2(J0*/3)CJJ  +2ll5«,  ^ 

"f  V«,., +  X  %,  (c3Tu-«w  ueo‘29o  ♦ 

-  C»J  -«o,[f  (flco«80-  8*in80)  ♦  co«80  -  7>s»m80)  -cj  <0 

♦  2  0b)  co*60+(  200 -8)sin$e]  +  C|SJ  ($  *0 )  *i  n0Q  ♦8cos0(J 

M*  0+  0  ( I  -CO *26J  ♦  -J-  { I  ♦  28  )  tin  20o]  cot2flyc*>(  -OTSCgjl-fljCgo, 

(1) 
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The  torsion*!  modal  squat  ion  from  Volume  I  of  this  report  it 


0* 


♦  C„q(  (i- Tj* -co*2S0) -CM  [?0+  -J-  ( 1  +  2i  ^  (i-coliS,)] 

-  II  q., g,e  C44i  #  *9,  II q. ,  * 9.n>, «vC>V«i 

-20&o+  &  Sin2®0  -  -y-  ( !♦  2&)(i-cos20o))  +  sm26()Cv5  j 

*•1 

c«,;  *s.,(  8°^  ')  'I9-}2*'  c*v  *■  xr 

*0 1  p  -co*J90l  ♦  -j-  (i  ♦  28 )  smzej  -  cot  Jfl0Cll4  ] 

-Il2C„,  q.  iq  co»90 -q  9*in90)+IIjc(,  q.  (q„  «.n9„ 

I  p  «•*!  •  9  9  v  i  j< 

-Q co*0o)-CM(  [e(£co*0o-  8tm^)  ♦  gl9co*0o-gy9»in0j 
-Cf3(  [  (5^/3^2>3d) cos«0-{5 -2/90)*m€^]  *C^[  {0  +  0)%\<\Oo 
♦Scot  90]  ♦Iq*j,CHS|  (/[  i0+0)%it>6o-  Scof^J  (2) 


The  flap  angle  equation  from  Volume  I  of  th's  report  is 


0  '  -Wfo^/lsA^d?-  MBrCQ(?0^g  ) -l9[S+0H  +  28)] 

-  C97[d0cosd0-  eos^eo-209os\r\Bo  -  {  \  -28  )sm0o] 

-  tonS^^I  (yfc'fc.oQ^  [Zq^  C.0  cosS, 
♦£^oCiSoS'nVSin*0C96  ] 

p  O 
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The  model  conatante  C,  through  C,,«  are  defined  in  Appendix  IV  of 
Volume  I.  They  are  reproduced  below  for  convenience 
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The  nondlmens  tonal  aerodynamic  twisting  moment  and  ths  non- 

d  las  ns  tonal  aerodynamic  shearing  fores  (5*)i,  appear  In  ths  above 
equations,  and  are  psrtlnent  to  ths  static  divergence  analysis.  Ths  ex* 
press  loos  for  these  quantities  are  given  as  Eqs.  ( lhl )  and  ( 1U6)  In  Voluae 
I  of  this  report.  They  are  reproduced  below  as  Eqs.  (5M  and  ($5). 


($»)«,  "  (*'*“♦  *  ctf*in+)  (5M 

♦yUtfH*  (CjCOS  <t>  ♦  Co  sm  6) 

-  tin  0  ( c,  sin  *  -  cdco*  *  )j|  ( 55 ) 

The  angle  •  is  given  by  Eqs.  Ub7)  through  ( 1*»9)  of  Voluae  1.  The 
aatheaatical  definition  of  •  Is  reproduced  below  as  Eq.  (56). 

®  “  075n  -  A,d  cos  f  -  B|$fin  qr  .  ton  8y  ^/9 +(Wf)r  ic 

♦  0,  (V  *  T  •  0  75)  ♦  0,  (  56) 

The  ai^tle  <p  is  Che  local  angle  between  the  plane  of  rotor  rotation  and 
the  local  resultant  velocity  vector,  as  shown  In  Figure  5  of  Volume  I. 

In  accordance  with  assumptions  3  and  7  of  this  voluae,  which  were  given  on 
page  5  ,  Eq.  (56)  becomes 


®  s  09  ~  0  tonBi  -  y ) 

Hote  that  8»  does  not  exist  for  the  nonart iculated  blade. 

In  accordance  with  assumption  2,  Eq.  (57)  becomes 

®  s  ^9,  00.  '  £  ton  ®»  '  ‘  J > 
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A m  in  Figure  5  of  Volume  I, 


«r  *  ®  ♦  ♦  (59) 

In  accordance  with  assumptions  k  and  8,  and  through  tha  usa  of  Eq  ($8) 
and  (59) 


(C|CO»*  ♦  Cgs.n*)  *  o(yg,Qg,  -  4»onl,  ♦  (50) 


is  obtained  for  conventional  flrw.  For  reverse  flow, 


(C|  cos  0  +  Cd  s»n^  )  *  -o(yg(qg,  -  $  long,  ♦  w  -  (51) 


is  obtained  In  the  same  way. 

When  Eq.  (60)  oi  (6l)  is  substituted  in  Eq.  (5*0.  the  expression  for  the 
nondinensional  flatwise  aerodynamic  shear  force  for  conventional  flow 
becomes 


(5.)z,  s  i(^5")c(lS fr)  0  *  r0,  Qg,  -  B  »°n  ( & ) 


and 

(5.).,  =  -  ±(<£)  c  (&)'  «  «  ’»,*>,  ■  *'«"»•  ♦  *  •  ♦  1  (63) 

for  reverse  flow. 

For  conventional  flow,  and  in  accordance  with  assumptions  1*,  5,  and  8, 

C"c/*  s  0  (61*) 

for  reverse  flow, 

Cmc/«  :  -  C  |  /  2  (65) 
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When  Iq.  (60)  and  (61*)  or  (61)  and  (6$)  are  uaed  in  Eq.  (55).  and  consid 
eration  la  given  to  aeaunptioni  1,  5.  and  7.  the  expression  for  the  non- 
diaenslonal  aerodynamic  twisting  mocae.it  becomes 


<"*>■.*  ±  (<£){*.« 1  (a fe)2  o  ( ’i, »«,- « '»»**♦) }  <<*> 

for  conventional  flow,  and 

OMs,  «  i  (-Sf)  {^*(iJt)2  f  •  *ton,»  ♦  *  -  #) 

-1*ca  C(^)*  o(rfliqe,-  /9»onl,^  *  -  #  )}  (67) 

for  reverse  flow. 

Eq.  (67)  can  be  written  in  the  following  form  by  collecting  terms. 

("Uh,*  £  {  (ife)2  (1  "  *•«*)  Co  (>M*i"0tan8*  +  v  ’$)}  (68) 

In  accordance  with  assumption  6. 

*  *  f°n$  *  (69) 

for  conventional  flow,  and 

ir  -  <p  2  ton(w-$)  s  -  (70) 


for  reverse  flow. 


The  expressi  on  for  UT/ilR  is  given  as  Eq.  (139)  of  Volume  I.  When  all 
dynamic  terns  are  eliminated  in  accordance  with  assumption  1,  the  follow¬ 
ing  expression  is  obtained. 

5r  s  m[('- •§■)  *in 'J'  ♦  Scos^]  +  f  +  v(i-f-)  * 


-  /?[We  cos 9  +(ve  +  y10lc/4)*infl] 

+  (wesin^  -  Vecos0)[-ft(cos^ -Ssin^)  -  Xs/3  -  y»o 5e/a co*©] 

-  0'[vesin0  +  Wfccosg]  ficot^f 

For  this  static  divergence  analysis,  using  assumptions  7  and  3, 

9 =  -/3  ton  g, 

9'  =  0 


(71) 


(72) 

(73) 


When  assumptions  2  and  7  are  used  on  Eq.  (71),  considering  Eqs.  (72)  and 
(73), 


Ml 

a* 


/isin^/  +  i  +•  f 


(7«0 


The  expression  for  UP/ilR  is  given  by  Eq.  (l^O)  of  Volume  I.  When  all 
dynamic  terms  are  eliminated  in  accordance  with  assumption  1,  the  follow¬ 
ing  expression  is  obtained. 

Sh  =  fcsO'f*)  -  pplcosy-  Ssint^)  -  /3  [( Ve  +  fo K/Jcos0  -^sinfi] 

+  (Wecos 9  +  V«' sing)  [-^(cosf  -  Ssin^)-  Xs/3  -  yioSC/4co*0] 

-  9'fJ-  cos^(Vecos^  -  ^esin^)  (75) 


When  assumptions  2  and  7,  and  Eqs.  (72)  and  (73)  are  employed  with  Eq.  v 7 5 X 
the  following  simplified  expression  for  Up/flR  results. 


Sr  =  x*  "  M/9cos^  -  yl0,c/4*e 


(76) 
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The  terns  containing  Viok/4  in  Eq.  (76)  are  small*  and  will  be  neg¬ 
lected  in  the  remaining  development. 

When  Eq.  (76)  is  used  with  Eqs.  (69)  and  (70)  and  assumption  6  in  Eqs. 

(62)  and  (63)*  the  following  expressions  result. 

The  expression  for  conventional  flow  is 

(50z,  s  £  (■flcr)c[(SSj2°(y«.^,  -/3fon8>) 

+  (Sr)  °(x»  -tipco%y  -/xWecos^)]  (77) 

The  expression  for  reverse  flow  is 

(5*)2#  *  [(sr)2  Q  (ty'fc,-  £,on8») 

“(5r)  °(x»-m£co»^ -/xWfCos*)]  (78) 

When  Eqs.  (69)*  (70),  (76),  and  assumption  6  are  used  in  Eqs.  (66)  and 
(68),  the  following  expressions  result: 

The  expression  for  conventional  flow  is 

0 m*)x*  s  i  (fi^)  C  0  Vto c/4  [(sr)2  ( y 9,  -  £ ,on  8» ) 

+  (tj*)(*s  -ft/Scos^r  -/xWe  cos*)]  (79) 

The  expression  for  reverse  flow  is 

(W*)I§  =  i  (?  _y»w)  [( Hr)2  ( q®<  ~  ^  »an  S, ) 

~ (sr)  C x*  ~  c°* ^  -  /x  We eo*^ )  ]  (so) 

The  algebraic  sign  of  Uj  as  given  by  Eq.  (7b)  determines  whether  a 
particular  blade  section  is  ?n  conventional  or  reverse  flow.  If  Ut>0 
the  blade  section  is  in  conventional  flow.  If  Uj<0  the  blade  section  is 
in  reverse  flow.  Therefore  a  single  expression  can  replace  Eqs.  (77)  and 
(78),  through  the  use  of  the  absolute  value  of  Ut  .  This  expression  is 

(5»)*#  =  i  ('^)^0[(flR)^'  tonSj) 

+  (X, -p/3co%+  ~  /xW^cos*)]  (81) 
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Similar  considerations  can  be  used  to  provide  a  single  expression  for 
,  also  valid  in  both  forward  and  reverse  flow,  from  Eqs.  (79) 
(80).  This  expression  is 


5  i  [ Tnfefi' ( "S 'on8i) 

♦(BDiw*  cos^ -/iWicos^)]  (8  2) 


Where  the  quantity  ▼sc  is  given  by 


7» 


CM 


-f(^ 


|Uth 


(83) 


The  flatwise  bending  modal  equation  was  presented  in  this  volume  as  Eq. 
(1).  This  equation  is  now  simplified  by  using  assumptions  1.  2,  3,  7* 
and  8  listed  on  page  5  .  The  resulting  equation  is 


0 


,T 

VS.)*,  df 


4 


dg((c»Tlt(  -CM(>1  -C..J  -  Cio,(Qtt) 


( 8U) 


The  Cbt,,.  modal  constant  accounts  for  spar  centroid  to  elastic  axis  dis¬ 
tance,  and  tip  supported  counterweight  effects.  This  modal  constant  is 
small  and  is  therefore  deleted  hereafter. 

The  torsional  modal  equation  was  presented  in  this  volume  as  Eq.  (2).  This 
equation  is  now  simplified  by  using  assumptions  1,  2,  3,  7,  and  8  on  page 
5  .  The  resulting  equation  is 

°  =  (^0  +  C*T<  +  ^so, 


-  C4«,  [-/3ton8i]  +  qW|  [-  Css,.,,  -  f CiS,t,  -  C#4,.,  ] 
-CM,  [f/3  +  g„]  -C„,[/9] 


(85) 
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The  CMli,  modal  eomtant  accounts  for  spar  centroid  to  elastic  axis 
distonce’ond  tip  supported  counterweight  effects.  This  constant  is  small, 
and  is  omitted  from  the  following  development.  The  term  Involving  C««, 
is  also  small,  and  is  deleted  in  the  following  development.  This  term  re¬ 
flects  the  steady  centrifugal  twisting  couple  due  to  blade  flapping  when 
pitch-flap  coupling  is  present. 

The  flap  angle  equation  was  presented  in  this  volume  as  Eq.  (3).  This 
equation  is  now  simplified  by  using  assumptions  1,  2,  3,  7  and  8  of  page 
5  .  The  resulting  equation  is 

0  8  (j i  <«.)..*")  *  +  -  Is* 


+  CpffitonSf  +  CNl£ton8| 


♦  QgOfe,)**  5jo  ,on®i  (86) 


The  terms  containing  CM  and  Ctr  are  small  compared  to  the  others,  and 
will  therefore  be  neglected  hereafter. 

When  the  articulated  blade  is  considered  in  accordance  with  assumption  3 
of  page  5  ,  the  following  equation  results  from  the  substitution  of  Eq. 
(82)  in  Eq.  (85)  : 


When  the  articulated  blade  is  again  considered  in  accordance  with  assump¬ 
tion  3  on  page  5  »  the  following  equation  results  from  the  substitution  of 


Eq.  (81)  in  Eq.  (86)t 
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(88) 


Mien  the  procedure  used  to  obtain  Eqs.  (8?)  and  (88)  is  repeated  for  the 
nonart iculated  blade,  Eq.  (82)  is  substituted  in  Eq.  (85),  and  Eq.  (8l)  is 
substituted  in  Eq.  (8M.  The  following  pair  of  equations  result. 


(89) 
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(90) 


In  the  above,  the  substitution 


has  been  made. 


(91) 


Eqs.  (87)  and  (88)  and  Eqs.  (89)  and  (90)  are  two  sets  of  tvo  linear 
equations  each.  These  equations  could  be  solved  for  qg,  and  /? 
from  Eqs.  (87)  and  (88)  for  articulated  blades,  and  for  q g,  and  qW|  , 
from  Eqs.  (89)  and  (90)  for  nonart iculated  blades,  if  all  the  other  par¬ 
ameters  were  given.  However,  the  objective  is  to  find  conditions  for 
boundaries  of  static  stability.  The  boundaries  of  static  stability  are 
those  sets  of  conditions  for  which  static  equilibrium  is  present  when  the 
right  hand  sides  of  Eq.  (87)  through  (90)  are  zero.  It  should  be  noted 
that  the  right  hand  sides  of  these  equations  represent  steady  torsional 
and  flapping  or  flatwise  bending  moments,  independent  of  the  deflections 
and  consistent  with  the  simplifying  assumptions  employed.  The  conditions 
for  a  static  stability  boundary  are  thus  those  conditions  such  that  in¬ 
definitely  large  static  deflections  can  occur  even  though  these  independ¬ 
ently  applied  steady  forces  are  zero. 


In  this  way,  the  following  pair  of  equations  was  obtained  for  articulated 

blades : 


0  = 


{■§&.(  V?«c° 51  loTi  d?  -CM^sjqe, 

+  £rr£  C o|uT|df  + 


7*c  Co  Urluil  dr 

iCf*  +  C*s^ 


(92) 
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0  = 


Ot|OtI  df  "  {ty,  )m  *o°8j  Sjo  ^  <J0( 

-  ,on®*^  ^Co0T|0rldF 

+  £Sl  cos* J  VCo|0tI  df  +  Mt  Fcg  I  ♦  Itj-  0 


(93) 


In  Eq.  (92)  and  (93)  t  the  tern  Ut  /SIR  has  been  replaced  by  0T  .  This 
approximation  is  in  considerate  of  assumptions  2  and  7  on  page  5  . 

In  a  similar  manner,  the  following  pair  of  equations  was  obtained  for  non- 
articulated  blades: 

0  =  •£  ^*c  Co  °TiaTi  df  ~ c* a®1}  % 

"  /*  C0,*J  y#'  WoW  dr  +  iC»>  +  c"}  «w,  ( 91* ) 


(95) 


In  Eq.  (92)  through  (95). 


-  ~  UT 

ut  *  SR"  (96) 

uT  =  9  ♦  f  +  fxiinf  (Q7) 
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z 


53  o 

T. 


GJo/mol  AR*)1 

(98) 

Tff,T*  df 

(9f>) 

A 

1  m  7  df 

(100) 

In  farther  developments,  the  following  more  convenient  notation  will  be  used 
for  the  distance  from  the  elastic  axis  to  the  center  of  gravity t 


(101) 


In  the  case  of  an  articulated  rotor,  the  determinant  of  the  two  homogeneous 
equations  (92)  and  (93)  is  formed  and  evaluated  for  many  values  of  advance 
ratio  ,  with  all  other  input  being  held  constant.  The  value  of  /x 
for  which  this  determinant  is  zero  is  found  within  the  desired  degree  of 
approximation  by  a  computer  trial-and-error  procedure.  This  value  of  /x 
is  considered  to  be  the  value  at  which  torsional  divergence  occurs.  The 
process  is  similar  with  a  nonarticulated  rotor  except  that  Eq.  (91*)  and  (95) 
are  used. 
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TORSIONAL  DIVERGENCE  DESIGN  CHARTS 


DISCUSSION 


The  limits  on  rotor  operating  conditions  imposed  by  blade  torsional  diver¬ 
gence  have  been  studied  extensively  by  means  of  the  relatively  simple  method 
described  above. 

A  fairly  large  number  of  parameters  have  been  selected  for  this  study,  with 
various  combinations  being  chosen  to  be  of  maximum  futui  e  use  in  the  pre¬ 
liminary  evaluation  of  new  blade  designs.  Design  charts  have  been  generated 
to  facilitate  such  evaluations  without  performing  new  calculations. 

The  various  parameters  criginaa.ly  chosen  to  be  varied  were  aa  follows: 

1.  Blade  inertia  ratio,  m0  k/  /p C0*oR*- 

2.  Control  system  spring  ratio,  KrLB/GJo* 

3.  Blade  frequency  parameter,  GJo/nr»0h/(HR)* 

L.  Blede  chord  distribution. 

5*  Blade  torsional  stiffness  distribution. 

6.  Blade  pitching  inertia  distribution. 

7.  Blade  average  lift-curve  slope. 

8.  Aerodynamic  root  cutout. 

9.  Blade  center-of-gravity  dictribution. 

In  addition  to  the  above,  the  effect  of  blade  lift-curve  slope  spanwlse 
variation  caused  by  Mach  number  effects  was  investigated.  The  effects  of 
pitch- flap  coupling  were  evaluated,  and  calculations  were  carried  out  for 
nonart icul at ed  blade  designs. 

Blade  designs  with  uniform  chord,  mass  and  stiffness  properties  were  used 
in  order  to  obtain  design  charts  for  the  effects  on  torsional  divergence 
boundaries  of  control  system  flexibility,  aerodynamic  root  cutout,  and 
center-of-gravity  position.  The  independent  effects  of  the  ether  parameters 
mentioned  above  were  studied  with  additional  calculations.  The  various 
charts  and  pertinent  assumptions  used  in  their  preparation  are  described  in 
this  subsection.  The  significance  of  the  simple  energy  analysis  results 
and  the  permissible  analytical  assumptions  and  preliminary  design  procedures 
are  given  later  under  the  appropriate  subsections ,  after  the  Normal  Mode 
Transient  Analysis  results  have  been  presented. 
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The  inertia  distribution  is  included  as  a  parameter  in  this  study,  even 
though  the  static  stability  calculation  does  not,  strictly  speaking,  con¬ 
sider  inertial  effects  as  such.  The  natural  vibration  mode  is  used  as  an 
assumed  divergence  mode;  thus.  Inertial  effects  on  the  divergence  calcula¬ 
tion  are  actually  a  reflection  of  changes  in  mode  shape.  Certain  unpub¬ 
lished  preliminary  results  have  indicated  that  correct  trends  are  predicted 
by  using  the  i  itural  mode  shape  in  this  manner. 

It  is  Instructive  to  consider  the  differential  equation  of  torsional  static 
deflection  of  a  rotor  blade.  This  differential  equation  expresses  the 
equilibrium  of  static  elastic  and  aerodynamic  torsional  couples  on  an  in¬ 
finitesimal  spanvise  increment: 


^(^3?)  5  "QCJC*  s  -qosC0  ( 10?) 

It  is  assumed  that  the  blade  is  uniform,  with  properties  of  the  reference 
75%  radius.  Then  the  differential  equation  is 


GJ  fjf  2  -qCxC0#0  *  -qo0«eC0  (103) 

The  symmetrical  blade  is  assumed  to  be  operating  at  zero  angle  of  attack, 
except  for  the  angle  caused  by  the  torsional  deflection. 

The  dynamic  pressure  q  due  to  velocity  components  normal  to  the  radius 
is  given  by  the  following: 


q  =  £  (Xlr  ♦  Vsin*)*  *  £  (i  ♦ 


:  id.) 


Nondimens ionali zing  and  using  Eq.  ( 10U )  in  Eq.  (103)  gives 

lfl$f  s  -f  (105) 


For  the  case  where  270°,  and  with  fx  >  I  ,  reverse  flow  exists  over  the 
entire  span,  and  e0/C0  equals  0.5  for  the  normal  25%  chord  center-of- 
gravity  and  elastic  axis  positions. 

Thus,  the  differential  equation  becomes 

$  5  f'R,(.-M>*g 


dot) 
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The  nondimensional  parameter  S„ 


Is  defined  by 


S. 
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Rigorous  s.  tisfaction  of  the  above  differential  equation,  assuming  fx  >  I  , 
is  provide'  by  the  following  infinite  series: 


0  =  P, 


^[(25^ 


Jl^L 


(4n)(4n-l)(4n-4)(4n-5)(4n-8)(4n-t)  (4)0) 
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P'  f“1l(2S#)i  (4«  *  I)  (4n)(4A  -3)(4< 


3)(4n-4)(4n  *7)(4n-t)  (SK4), 


(109) 


The  series  given  as  Eq.  (109)  converges  fairly  rapidly,  so  only  a  few  terms 
need  be  used  for  numerical  evaluations. 

The  appropriate  boundary  conditions  for  the  blade  fixed  at  the  root  are 


e  -  o 

when  i :  0 

( 110) 

o 

when  s  s  | 

(111) 

C-e  of  these  with  Eq.  (109)  provides  two  homogeneous  equations.  Setting 
the  determinant  equal  to  zero  provides  a  characteristic  equation  in  i/S» 
if  ^  is  given.  The  lowest  value  of  I/S*  which  makes  the  determinant 
zero  is  the  critical  value  of  interest.  When  this  has  been  found,  the 
value  of  pt  in  terms  of  p0  can  be  found  by  using  one  of  the  boundary 
conditions.  The  above  procedure  was  carried  out  for  p  ~  i  with  tho  results 
shown  in  Figure  1  and  with  a  solution  for  critical  S*  ■  .031.  Also  shown 
for  comparison  is  the  uniform  blade  m.tural  vibration  mode. 

The  digital  computer  was  used  to  predict  a  divergence  advance  ratio,  using 
Eq.  (92)  and  (93)  as  described  previously.  The  first  natural  torsional 
vibration  mode  was  used  in  these  calculations.  With  a  value  of  S*  *  03i  , 

the  predicted  value  of  critical  advance  ratio  was  /u  ■  1.0J,  compared  with 
fX  ■  1.0  for  the  differential  equation  solution.  This  demonstrated  that 
the  use  of  a  natural  vibration  mode  as  a  divergence  node  ic  a  good  simplify¬ 
ing  assumption,  even  for  advance  ratios  near  unity. 


The  ao«t  important  conclusion  to  be  drawn  from  the  above,  however,  is  that 
the  advance  ratio  ft,  and  the  torsional  stiffness  coefficient  S»  #  are 
the  basic  paraseters  for  the  study  of  retreating  uniform  blade  static  tor¬ 
sional  divergence.  The  same  parameters  also  arise  from  inspection  of  a 
single-degree-of-freedoo  energy  analysis,  where  S»  apj ears  as  a  parameter 
along  with  the  nondimens ional  integrals  of  strain  enerry  and  virtual  work. 
The  parameter  S*  is  also  the  product  of  blade  inert!*-  ratio  and  frequency 
parameter,  divided  by  average  lift-curve  slope.  These  parameters  arc  among 
those  listed  at  the  beginning  of  this  section.  Variations  in  these  par¬ 
ameters  are  included  as  variations  in  the  more  basic  S»  parameter.  The 
numerical  value  of  S«  for  a  typical  rotor  blade  operating  at  design  tip 
speed  is  approximately  (X0l»0. 

The  parameter  S«  has  been  developed  with  a  $0%  chord  distance  between  the 
aerodynamic  center  of  pressure  and  the  center  of  gravity,  which  is  coin¬ 
cident  with  the  elastic  axis.  This  is  appropriate  to  the  usual  retreating 

blade  case. 

The  usefulness  of  attempting  the  definition  of  a  parameter  even  more  gen¬ 
erally  applicable  than  Sa  has  been  examined.  Specifically,  Eq.  (106)  can 
be  rewritten,  retaining  the  nondlmens ional  distance  between  the  aerodynamic 
center  of  pressure  and  the  center  of  gravity,  which  is  coincident  with  the 
elastic  axis. 


g  •  M.  ♦*»»♦>*  8  (U2! 

It  is  evident  from  consideration  of  the  right  hand  side  of  Eq.  (102)  that 
divergence  will  first  appear,  as  aircraft  forward  speed  increases,  at  that 
azimuth  angle  where  the  coefficient  of  0  reaches  a  maximum  absolute  value. 
Prom  Eq.  (10M  it  can  be  seen  that  two  particular  azimuth  angles  are  of 
primary  Interest;  these  are  *  ■  90°  and  +  »  270°.  Note  that  the  value 
of  So/C«  is,  in  general,  different  at  these  two  azimuth  angles.  For 
*  ■  90°,  for  example 


(113) 


while  for  +  ■  270°,  and  fi  >  l  , 

-  50  ♦  Vc*  / iij ) 

where  ?ee  1b  cons  ent  along  the  blade  radius.  When  only  the  two 
specific  azimuth  angles  ♦  ■  90°  and  ■  270°  are  considered,  Eq.  (112) 

has  two  possible  forms.  These  are 
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for  \p  =  270°.  It  night  be  concluded  from  the  above  that  the  coefficient 
of  (*+^)2e  and  (x  > l)zQ  in  the  above  equations  is  the  most  general  par¬ 
ameter.  This  would  oe  true  if  only  advance  ratios  greater  than  unity  were 
of  interest.  For  advance  ratios  ft  <  1,  the  above  simple  concepts  are 

not  true  even  for  a  uniform  blade,  since  e0/C0  will  vary  radically  between 
the  inboard  portions  of  the  blade  where  reverse  flow  is  present,  and  the 
outboard  portions  where  conventional  flow  is  pres-nt.  An  effective  value 
of  e0/C0  might  be  defined,  but  this  itself  would  be  a  function  of  advance 
ratio.  Therefore  it  is  advantageous  to  accept  the  Sn  parameter  as  the 
basic  parameter,  and  consider  variations  in  e0/C0  through  the  generation 
of  both  advancing  and  retreating  blade  stability  boundaries,  and  through 
the  treatment  of  Yca  as  an  independent  parameter. 

The  energy  analysis  described  during  the  development  of  Eq.  (92)  through 
(95)  introduces  a  flapping  or  bending  degree  of  freedom  to  the  torsional 
divergence  analysis.  This  accounts  for  the  fact  that  blade  static  torsional 
stability  always  involves  some  blade  flapping  due  to  lift  caused  by  tor¬ 
sional  deflection.  The  resistance  to  flapping  of  conventional  rotor  blades 
is  provided  mainly  by  a  centrifugal  force  component  which  acts  at  the  blade 
chordwise  center  of  mass. 

If  the  center  of  mass  does  not  coincide  with  the  aerodynamic  center  the 
resulting  couple  causes  potentially  divergent  torsional  deflections.  At 
azimuth  angles  of  270°  and  90°,  the  results  for  articulated  blades  with  no 
pitch- flap  coupling  are  practically  identical  with  single-degree-of-freedom 
analyses  if  the  torsional  couple  caused  by  the  distance  between  the  aero¬ 
dynamic  center  and  the  center  of  mass  is  considered. 

The  uniform  blade  design  charts,  Figures  2,  3,  and  1,  were  produced  for  a 
completely  uniform  blade  with  a  chcrd-to-radius  ratio  of  0.050  and  negligible 
flapping  hinge  offset.  The  value  of  mass  ratio  Mr  was  11.9.  Calculations 
were  carried  out  for  practical  variations  in  this  parameter,  and  divergence 
results  were  identical  within  the  reading  and  data  input  accuracy.  The 
ratio  of  the  blade  torsional  radius  of  gyration  to  rotofc  radius  was  0.012. 

The  elastic  axis  was  assumed  to  lie  along  the  25/5  chord  line.  The  lift- 
curve  slope  was  assumed  to  be  constant  along  the  blade  span,  with  no  tip 
loss  being  considered. 
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The  boundaries  were  generated  by  finding  the  critical  advance  ratio  for  10 
different  values  of  nondinsnsional  stiffness  coefficient  Sr  for  each  of 
the  otherwise  identical  blade  configurations. 

Torsional  elastic  stiffness  was  changed  to  vary  the  nondimens ional  stiff¬ 
ness.  Identical  values  of  torsional  stiffness  coefficient  Sr  were  gen¬ 
erated  by  changing  rotational  speed  with  elastic  stiffness  held  constant, 
and  calculations  were  carried  out  to  verify  that  identical  values  of  crit¬ 
ical  advance  ratio  would  be  obtained.  Sms 11  variations  in  torsional  mode 
shape  due  to  different  rotational  speeds  were  not  found  to  be  significant. 
The  effects  of  rotation  on  blade  torsional  stiffness  were  also  unimportant. 

The  first  set  of  design  charts,  Figure  2,  presents  plots  of  the  critical 
advance  ~atio/x  versus  torsional  stiffness  coefficients  Sr  for  various 
values  of  aerodynamic  root  cutout  Xqa  •  Each  plot  is  for  constant  values 
of  control  stiffness  parameter,  l/K  ,  and  chordwise  center-of-gravity 
offset,  Yco/C  .  These  plots  are  for  critical  values  of  occurring  at 

azimuth  singles  of  either  270°  or  90°.  The  90°  azimuth  torsional  divergence 
occurs  only  for  the  aft  center-of-gravity  positions  (negative  YCo/C  ), 
and  is  represerted  by  the  branch  on  the  left  side  of  the  break  in  these 
curves.  The  divergence  boundaries  for  the  90°  azimuth  are  shown  only  as 
they  are  more  critical  than  the  270°  azimuth  boundaries,  and  vice  versa. 

This  set  of  curves  will  be  useful  in  the  evaluation  of  specific  blade 
designs. 

The  second  set  of  design  charts.  Figure  3,  presents  cross  plots  of  the 
first  set  at  various  constant  values  of  p  and  Yc#/C  .  The  values  of  ^ 
presented  are  1.2,  l.U,  1.6,  1.8,  and  2.0,  which  cover  the  range  expected 
to  be  of  greatest  future  interest  for  torsional  divergence.  Both  positive 
and  negative  values  of  Ycr/C  are  included.  The  values  of  YCq/C  are  .10, 

0,  -.02,  -.06,  and  -.10.  The  reciprocal  of  the  control  stiffness  ratio  K 
is  plotted  against  Sr  at  constant  values  of  Xqa  •  This  set  of  curves 
can  be  used  in  trade-off  studies  concerning  the  relative  benefits  of  changes 
in  Sr  and  K.  Since  Figure  3  is  composed  of  direct  cross  plots  of  Figure 
2,  it  also  reflects  instabilities  at  azimuth  angles  of  either  90°  or  270°, 
depending  on  which  is  most  critical. 

The  third  set  of  design  charts,  Figure  l*,  also  presents  cross  plots  of  the 
first  set,  but  with  Yc«/C  plotted  against  Sr  at  various  constant  values 
of  Xqa  •  The  values  of  ft  and  l/K  remain  constant  for  each  of  these 
plots.  This  set  of  curves  can  be  used  in  trade-off  studies  concerning  the 
relative  benefits  of  changes  in  Sr  and  YCo/C  .  As  with  Figure  3,  the 
boundaries  are  shown  for  azimuth  positions  of  either  90°  or  270°,  depending 
on  which  i3  most  critical.  The  azimuth  angle  to  which  a  particular  point 
on  any  of  the  boundaries  pertains  is  easily  determined  by  reference  to 
Figure  2. 

The  independent  effects  of  certain  other  parameters  on  the  static  divergence 
stability  boundaries  are  shown  in  Figures  5  through  lU.  The  value  of  Sr 
in  these  cases  was  determined  by  using  the  values  of  GJo  ,  0o  ,  and  Co 
existing  at  the  75%  radius.  The  shapes  and  frequencies  of  each  of  the  tor¬ 
sional  modes  are  appropriate  to  each  configuration.  The  rigid  flapping 
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mode  was  used  for  deflections  out  of  the  plane  of  rotation  for  the  articu¬ 
lated  rotor.  The  first  flatwise  bending  mode  appropriate  to  the  given 
values  of  a  parameter  m/^OoCcT,  the  ratio  of  torsional  to  flapwise  bending 
stiffness  EIj/GJo  »  and  the  torsional  stiffness  coefficient  Sr  was  used 
in  the  case  of  the  nonart iculated  blades.  In  eacn  case  where  blade  cross- 
sectional  properties  were  varied  separately  for  this  set  of  charts,  they 
were  made  equal  to  the  uniform  blade  properties  at  the  75 ?  rotor  radius. 
Information  as  to  how  the  various  parameters  were  varied  is  contained  on 
the  captions  and  legends  for  each  figure. 

The  effect  on  retreating  blade  torsional  divergence  of  a  linearly  varying 
blade  chord  is  shown  in  Figure  5*  for  different  values  of  tip-chord-to-root- 
chord  ratio.  The  chord-to-radius  ratio  at  the  75?  radius  remained  at  the 
value  used  for  the  uniform  blade,  and  Sr  is  also  based  on  this  value. 

All  other  blade  properties  remained  as  for  the  uniform  blades. 

The  effect  of  a  structural  flexibility  cutout,  X0t  on  retreating  blade 
torsional  divergence  Is  shown  in  Figure  6.  The  structural  flexibility 
cutout  refers  to  a  torsionally  rigid  inboard  section.  Thus,  when  X0#  =  .10 
the  inboard  10?  radius  of  the  blade  is  considered  to  be  torsionally  rigid. 
The  remainder  of  the  blade  is  identical  to  the  previously  considered  uniform 
blades.  The  torsional  deflection  modes  used  in  these  calculations  are 
natural  vibration  modes  appropriate  to  these  blades  with  torsionally  rigid 
inboard  sections. 

The  effect  on  retreating  blade  torsional  divergence  of  a  linearly  varying 
blade  torsional  stiffness  is  shown  in  Figure  7  for  different  values  of  tip- 
stiffness-to-root-stiffness  ratio.  The  nondimensional  stiffness  coefficient 
SR  is  defined  at  the  three-quarter  radius.  The  blade  torsional  deflection 
modes  used  in  these  calculations  are  the  appropriate  natural  vibration 
modes.  The  significance  of  the  small  apparent  effect  of  tapering  stiffness 
will  be  discussed  later. 

The  effect  of  inertial  cutout  on  retreating  blade  torsional  divergence  is 
shown  in  Figure  8.  Inertial  cutout  refers  to  a  uniform  blade,  with  an  in¬ 
board  section  of  negligible  torsional  inertia.  This  has  an  effect  on  tor¬ 
sional  divergence  through  the  use  of  the  natural  vibration  mode  as  an 
assumed  divergence  mode.  Thus,  comparison  with  the  completely  uniform 
blades  reflects  only  the  effect  of  the  change  in  mode  shape  on  the  static 
stability  calculation.  Except  for  the  cutout,  the  blades  were  identical  to 
the  previously  considered  uniform  blades. 

The  effect  of  linearly  varying  inertia  on  retreating  blade  torsional  diver¬ 
gence  is  shown  in  Figure  9-  The  inertia  was  made  equal  to  that  for  the 
previously  considered  uniform  blades  at  the  75?  radius,  with  different 
ratios  of  root  inertia  to  tip  inertia.  Except  for  the  variation  in  inertia, 
the  blades  were  identical  to  the  previously  considered  uniform  blades.  It 
should  be  noted  that  the  effect  on  the  torsional  divergence  static  stability 
stems  only  from  the  alterations  in  the  torsional  natural  mode  shape. 
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Torsional  divergence  cases  were  considered  for  the  retreating  blade  with 
apanvise  variations  in  the  chordwise  center-of-gravity  offset  position. 
These  included  offsets  of  Yc«/C  ■  .10,  extending  from  inboard  fractional 

radii  of  0.25  and  0.L5  to  the  tip,  and  linearly  varying  offset  distances, 
with  Yc*/C  equal  to  0.10  at  the  15%  radius.  Linearly  varying  offsets, 
such  that  the  ratio  of  root-to-tip  offsets  was  two  and  three,  were  chosen 
and  also  used  in  separate  sets  of  divergence  calculations.  The  results  for 
all  of  these  calculations  were  within  1  percent  of  the  values  for  the 
uniform  blade;  therefore,  no  separate  charts  were  made. 

The  next  charts,  Figures  10  and  11,  show  the  effect,  of  the  distribution  of 
Clauert's  lift-curve  slope  correction  due  to  Mach  number.  The  effects  on 
the  advancing,  as  well  as  on  the  retreating,  blades  are  shown,  "he  par¬ 
ameter  Sc,  defined  in  the  table  of  symbols,  is  als  equal  to  the  parameter 
Sr  ,  evaluated  with  a  lift-curve  slope  of  2  7 r  ,  and  multiplied  by  the 
square  of  the  ratio  of  rotational  tip  speed  ilR  to  the  speed  of  sound. 

The  value  of  Sr  in  Figures  10  and  11,  however,  refers  to  the  value  at 
the  15%  radius  with  the  lift-curve  slope  correction  applied.  The  zero 
value  of  Sc  refers  to  the  earlier  uniform  blade  cases  without  Glauert's 
correction.  A  0.01^3  value  of  Sc  refers  to  conventional  blade  stiffness 
coefficient  and  rotational  speed  under  the  standard  conditions,  as  typified 
by  the  CH-3C  rotor.  Iv  should  be  noted  that  the  curves  of  constant  Sc  in 
Figure  10  are  very  nearly  curves  of  constant  advancing  blade  Mach  number  at 
the  15%  radius.  The  advancing  blade  15%  radius  Mach  number  along  the  curve 
Sc  *  .0286  ranges  between  0.68  and  0.71. 

The  next  charts,  Figures  12  and  13,  show  the  comparative  divergence  bound¬ 
aries  for  articulated  and  nonarticulated  blades,  with  the  nonart iculated 
uniform  blades  rigidly  built  into  the  hub  art  the  center  of  rotation.  Ad¬ 
vancing  as  well  as  retreating  blade  divergence  is  considered .  Note  that 
the  nondimens ional  stiffness  concept  remains  valid  for  the  nonarticulated 
blades,  so  long  as  the  ratio  EIf/oo0C0*  /  m  GJo  remains  constant.  This 

point  will  be  discussed  later.  The  value  of  EIf/GJ0  of  Q  5^8  refers  to 
conventional  practice  as  typified  by  the  CH-3C  rotor. 

The  effect  of  pitch-flap  coupling  as  predicted  by  the  simple  divergence 
analysis  is  shown  in  Figure  lL.  The  effects  on  the  retreating  blade  only 
are  shown,  since  pitch-flap  coupling  will  greatly  increase  advancing  blade 
divergence  advance  ratio  and  decrease  it  for  the  retreating  blade.  The 
effects  o**  pitch-flap  coupling  were  investigated  for  the  articulated  blade 
only,  since  the  provision  of  an  appreciable  amount  of  this  type  of  coupling 
is  mechanically  difficult  with  nonarticulated  blades. 


30 


10RSI3NAL  DIVERGENCE  CALCULATIONS  WITH  THE  EXTENDED 
NORMAL  MODE  TRANSIENT  ANALYSIS 


DISCUSSION  OF  PURPOSE  AND  METHODS 


A  number  of  solutions  of  the  flexible  blade  equations  of  motion  for  tor¬ 
sional  divergence  conditions  were  accomplished,  by  use  of  the  extended 
Normal  Mode  Transient  Analysis.  The  solutions  demonstrate  the  effect  on 
torsional  divergence  '  f  blade  flatwise  and  edgewise  flexibility,  in  con¬ 
junction  with  consideration  of  all  significant  inertial  terms.  This  was 
done  for  coincident  and  noncoincident  section  elastic  and  center-of-gravity 
axes  and  for  articului.ed  and  nonarticulated  blades.  In  addition  to  this, 
cases  were  run  at  red  ced  air  density,  to  provide  data  at  higher  values  of 
torsional  stiffness  <:•  'efficient  SR  .  Control  system  flexibility  was  also 
investigated  for  the  articulated  blades. 

The  effects  of  pitch-flap  coupling  on  the  torsional  divergence  of  an  artic¬ 
ulated  blade  were  calculated,  and  the  effect  of  a  representative  system  of 
gyro  feedback  on  a  nonarticulated  blade  was  determined. 

The  solution  of  the  blade  equations  of  motion  was  generally  conducted  by 
starting  a  typical  blade,  whose  properties  are  given  in  Table  I  and  II  at 
an  azimuth  angle  of  zero.  Collective  pitch  and  inflow  were  also  made  equal 
to  zero. 

The  elastic  modes  considered  for  the  articulated  blades  were  the  first 
through  the  third  flatwise  bending  modes,  the  first  edgewise  bending  mode, 
and  the  first  torsional  mode.  These  elastic  modes  were  in  addition  to  the 
rigid  blade  flapping  and  lagging  modes  for  the  articulated  blades.  For  the 
nonarticulated  blades,  the  first  through  the  fourth  flatwise  bending  modes, 
the  first  and  second  edgewise  bending  modes,  and  the  first  torsional  mode 
were  considered.  The  solution  for  blade  motion  was  then  obtained  for  one 
revolution  of  the  rotor  by  using  the  extended  Normal  Mode  Transient  Analysis. 
The  values  of  aircraft  forward  speed  and  rotor  rotational  speed  at  the  var¬ 
ious  values  of  advance  ratio  were  chosen  by  setting  the  advancing  blade  tip 
Mach  number  equal  toO.  85.  The  aerodynamic  data  were  appropriate  to  an  NACA 
0012  airfoil,  with  Mach  number  and  stall  effects  included.  Rotor  shaft 
angle  of  attack  was  set  constant  at  zero  degrees.  The  parameter  used  to 
evaluate  blade  response  was  one-half  peak-to-peak  stress  variation. 

It  will  be  noted  that  the  effect  of  blade  center-of-gravity  offset  was 
studied  at  a  constant  value  of  advance  ratio,  with  a  collective  pitch  @.75* 
of  2  degrees  and  an  inflow  X  of  0. 

This  alternate  method  was  necessary  to  demonstrate  the  existence  of  the 
practical  operating  limits  appropriate  to  the  torsional  flexibility  of  the 
advancing  blade.  The  excitation  of  the  unloaded  advancing  blade  was  found 
to  be  very  small  at  the  predicted  static  stability  boundary.  Extensive 
penetration  of  this  boundary  witnout  important  blade  response  could  occur 
for  the  completely  unloaded  rotor.  However,  the  addition  of  a  small  load 
resulted  in  a  torsional  stress  rise  similar  to  that  for  t  e  retreating  blade. 
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In  view  of  the  objective  of  correlating  the  simple  energy  analysis  for  tor¬ 
sional  divergence  with  practical  operating  limits,  the  advancing  blade 
cases  were  carried  out  with  the  moderate  loadings  generated  by  the  2  degree 
collective  pitch  setting.  One  possible  explanation  for  the  lack  of  exci¬ 
tation  on  the  advancing  blade  is  the  stabilizing  center-of-pressure  move¬ 
ment  at  high  subsonic  Mach  numbers.  This  would  tend  to  prevent  torsional 
divergence  until  a  relatively  large  section  of  the  rotor  disc,  including 
the  lower  Mach  number  areas,  is  unstable.  A  small  amount  of  cyclic  pitch 
was  used  to  eliminate  rotor  pitching  and  rolling  moment  for  the  nonartic- 
ulated  rotor.  The  method  of  solution  was  otherwise  identical  to  the  other 
Normal  Mode  Transient  Analysis  cases  for  this  section,  with  the  starting 
values  for  the  blade  differential  equations  of  motion  being  given  as  zero 
flapping,  lagging,  and  elastic  modal  displacements  and  velocities  at  an 
azimuth  angle  <j/  of  0  degrees . 

The  pitch-flap  coupling  effects  were  evaluated  with  the  unloaded  rotor  as 
described  previously.  Inclusion  of  this  type  of  coupling  in  the  calcula¬ 
tions  is  a  routine  provision  of  the  Normal  Mode  Transient  Analysis. 

The  study  of  the  effect  of  a  control  gyroscope  on  torsional  divergence  is 
much  more  involved.  The  scope  of  this  work  is  such  that  one  typical  gyro¬ 
scope  installation  could  be  included.  The  control  gyro  parameters  of  in¬ 
terest  are  given  in  the  descriptions  of  configurations  to  follow. 

At  this  point  in  the  discussion,  it  is  appropriate  to  outline  the  manner  in 
which  the  gyro  is  assumed  to  act.  Consistent  with  the  remainder  of  the 
investigation.no  shaft  or  fuselage  motion  is  included.  In  addition  to  the 
torsional  moment,  a  small  component  of  blade  flapping  moment  at  each  of  the 
blade  roots  is  reacted  by  the  pushrod,  because  of  a  small  built-in  angle 
between  the  blade  axis  and  the  feathering  axis.  This  effect  has  been  found 
to  be  small  for  this  torsional  divergence  study,  since  bending  moments  are 
initially  small  for  the  unloaded  blade.  The  blade  moment  about  the  feather¬ 
ing  bearing  is  provided  by  the  gyro  reaction  to  pushrod  force,  acting 
through  its  moment  arm.  The  gyro  motion,  and  the  resulting  pitch  position 
input  to  the  blades,  depends  on  the  pushrod  loads  applied  by  all  the  blades. 

In  order  to  assess  the  effect  of  the  gyroscope,  blade  response  calculations 
were  first  performed  with  the  gyroscope  assumed  as  being  fixed  to  the  shaft. 
The  loads  on  the  gyroscope  were  calculated,  and  then  its  response  to  these 
loads  was  calculated.  Finally,  the  gyro  response  was  used  to  define  new 
cyclic  pitch  control  positions  for  another  blade  response  calculation.  The 
change  in  blade  response  due  to  the  new  control  positions  was  then  noted. 

The  process  could  be  repeated  to  improve  accuracy,  although  one  cycle  was 
Judged  to  be  sufficient  for  the  purposes  of  this  study. 

CONFIGURATIONS  AND  FLIGHT  CONDITIONS  STUDIED 

The  pertinent  characteristics  of  the  basic  rotor  system  used  for  torsional 
divergence  studies  with  the  extended  Normal  Mode  Transient  Analysis  are 
shown  in  Table  I. 
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The  chord,  stiffness,  mass,  blade  center- of-gravity  chordvise  position,  and 
section  properties  of  the  blade  are  constant,  except  for  the  aerodynamic 
and  structural  cutouts.  The  elastic  axis  was  assumed  to  lie  along  the  25Jf 
chord.  The  blades  had  an  NACA  0012  airfoil  section.  The  nonart iculated 
blade  was  identical  to  the  art iculated  blade  except  for  the  removal  of  the 
flap-lag  hinge.  Control  gyro  parameters  of  interest  are  a  lead  angle  of 
1.5°  between  blade  axis  and  feathering  axis,  and  a  gyro  polar  moment  of 
inertia  Cl 01  times  that  of  the  rotor.  The  pushrod  was  assumed  to  lead  the 
blade  feathering  axis  by  an  azimuth  angle  of  U 5° •  No  mechanical  damping 
was  applied  to  gyro  motion. 

The  various  configurations  used  for  this  part  of  the  torsional  divergence 
study  with  the  extended  normal  mode  analysis  are  shown  in  Table  II. 

The  various  flight  conditions  used  for  this  part  of  the  study  are  shown  in 
Table  III.  The  advancing  blade  tip  Mach  number  remained  constant  at  M,it0 
=  .85. 

The  combinations  of  configurations  and  flight  conditions  are  shown  in 
Table  IV. 

RESULTS  OF  NORMAL  MODE  TRANSIENT  ANALYSIS  CALCULATIONS 


Figure  15  shows  a  three-revolution  time  history  of  blade  first  torsion 
modal  response,  after  being  started  with  zero  collective  pitch,  cyclic 
pitch,  flapping,  lagging  and  elastic  modal  displacements  and  velocities,  at 
an  azimuth  of  zero  (directly  aft).  This  figure  is  for  the  basic  articulated 
rotor  configuration,  operating  at  an  advance  ratio  of  1.6.  Note  that  the 
retreating  blade  torsional  response  is  more  severe  during  the  first  revo¬ 
lution  than  during  the  following  two,  which  Justifies  the  use  of  a  single 
revolution  to  evaluate  peak  torsional  divergence  response. 

Figure  16  shows  a  similar  two-revolution  time  history  of  blade  first  tor¬ 
sion  modal  response,  with  the  15^  chord  aft  center-of-gravity  position, 
which  is  the  articulated  blade  configuration  6.  The  response  is  calculated 
with  2  degrees  of  collective  pitch  present  to  supply  additional  loading,  as 
discussed  earlier.  Note  that  in  this  case  the  first  and  second  revolution 
maximum  response  is  virtually  identical,  again  supporting  the  use  of  the 
first  revolution  to  evaluate  torsional  response  for  the  lightly  loaded 
blade. 

Figure  17  shows  one-half  peak-to-peak  stresses  calculated  during  the  first 
revolution  after  starting  the  blade  at  \f/  =  0  ,  with  no  initial  loadings, 
at  various  values  of  advance  ratio.  This  is  for  the  basic  articulated 
rotor.  Also  shown  is  the  stability  boundary  predicted  by  using  Figure  6. 
F'.gure  6  was  used  by  finding  the  intersection  of  the  X0s  =  .06  stability 
boundary  and  a  plot  of  versus  S*  appropriate  to  standard  sea  level 
atmospheric  conditions  and  the  advancing  blade  tip  Mach  number  ofO.85.  The 
method  of  obtaining  this  relationship  will  be  explained  later,  under  the 
general  discussion  of  preliminary  design  methods. 
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Figure  Id  is  similar  to  Figure  17,  with  calculations  carried  out  at  a  value 
of  air  density  p  at  2/3  of  that  for  standard  sea  level  conditions.  This 
corresponds  to  an  altitude  of  approximately  13,000  feet. 

Figure  19  shows  the  Normal  Mode  Transient  Analysis  one-half-peak-to-peak 
stress  results  for  the  basic  rotor  with  control  system  flexibility.  The 
effects  of  control  system  flexibility  on  the  static  torsional  divergence  of 
the  blade  with  structural  cutout  was  not  determined  with  the  simple  anal¬ 
ysis,  since  structural  cutout  was  a  parameter  whose  independent  effect  was 
determined.  In  order  to  estimate  the  static  divergence  boundary  for  the 
combined  effects  of  structural  cutout  and  control  system  flexibility,  the 
independent  effects  of  each  were,  in  effect,  superimposed.  This  procedure 
is  not  rigorous  in  a  strict  mathematical  sense,  but  is  Justified  by  the 
approximations  inherent  in  the  static  torsional  divergence  analysis.  The 
procedure  will  be  described  mathematically  in  the  following  development. 
This  particular  discussion  is  limited  to  the  case  of  the  blade  with  uniform 
spanwise  properties  except  for  cutout,  in  order  to  present  the  procedure  in 
as  simple  a  manner  as  possible. 

It  can  be  observed  that,  in  principle,  the  advance  ratio  for  stati  :  diver¬ 
gence  can  be  expressed  as  a  function  of  all  the  variables  considered. 


M  =  *  (S*  .  I/K  ,  X0,  ,X0A  ,  Yc.) 


(117) 


If  p.  is  known  for  a  particular  set  of  values  of  S„,  l/K,  X,* ,  XM  , 

and  a  Taylor's  series  may  be  written  about  that  value.  The  elements 

of  that  set  are  identified  below  by  the  "a"  subscript. 


+ 


(118) 


When  values  of  the  variables  considered  above  remain  suitably  close  to 
Sr,,  l/Ka  ,  Xqj^,  Xo^ ,  and  Yce  ,  the  specific  terms  appearing  in  Eq. 

(118)  furnish  a  satisfactory  approximation.  The  derivatives  in  Eq.  (118) 
are  evaluated  at  S*a,  l/K0,  Xos#,  X0A# ,  and  Vcg0 

The  particular  case  at  hand  will  now  be  treated  by  using  Eq.  (ll8).  The 
values  of  the  variables  with  the  "a"  subscript  are  chosen  to  be  consistent 
with  the  available  data.  They  are 


pj 

=  1.39 

(119) 

S.o 

=  .107 

(120) 

3^ 


It  is  desired  to  find 


l/Kfl  *  0 

(121) 

Xoio  s  0 

(122) 

X<3^  =  0 

(123) 

^COq  =  0 

(122) 

/x  for  the  values  of  the  variables  given  as 


l/K 

=  2 

(125) 

X0g 

=  .06 

(126) 

Xo* 

=  .15 

(127) 

Yea 

=  0 

(128) 

Note  that  ?co 

is  equal  to  ?co 

,  which 

causes  the  last  term  in  Eq. 

(118)  to  vanish. 

The  values  of  ixQ 

and  S. 

were  found  from  Figure  2, 

at  the  intersection  of  the  appropriat 

;e  divergence  boundery  and  a  plot  of 

/x  versus  Sr 

calculated  for  the 

advancing 

blade  tip  \'ach  number  of 

O.85.  This  particular  relationship  will  be  defined  later  in  Eq.  (lUl).  In 
the  present  case,  Eq.  ( lUl )  provides  a  relationship  between  p  and  S* 
which  is  written  as  the  following  equation. 

S*  =  OI855(  I +/x)2  (129) 

The  values  of  the  derivatives  in  Eq.  (ll8)  are  found  from  Figures  2  and  6 
by  dividing  an  increment  in  fi  by  an  increment  in  the  variable.  For 
example , 

&H-  ~  &£. 

aft)  Aft)  (130) 

where  the  increments  A/z  and  A  <I/K)  are  taken  while  all  other  variables 
are  held  constant.  The  values  of  the  derivatives  needed  are 


(131) 


*1L 

d{i) 


(132) 
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(133) 


Je. 

a(Xo*) 


4f  =  -08 


(13M 


When  the  values  given  by  Eq.  (119)  through  (128)  and  by  Eq.  (.131)  and  (13M 
are  used  in  Eq.  (ll8),  the  following  equation  is  obtained: 


M  S  1.30  ♦  3.2  (S*  -  .107)  =  98  +  3.2S* 

(135) 

When  Eq.  (135) 
found.  These 

is  used  in  Eq.  (129)  the  final  values  of  p 
values  are 

and  S  *  are 

fi  =  1.30 

(136) 

S*  =  .10 

(137) 

The  effects  of  pitch-flap  coupling  are  reflected  in  Figure  20  for  the  re¬ 
treating  blade.  The  combined  effects  of  pitch- flap  coupling  and  structural 
root  cutout  on  the  static  torsional  stability  boundary  were  found  by  using 
Figures  6  and  ll*  with  a  method  analogous  to  that  demonstrated  by  Eq.  (117) 
through  (137). 

The  effects  of  aft  center  of  gravity  on  the  articulated  blade  are  shown  in 
Figure  21  for  an  advance  ratio  of  Q60.  The  rotor  was  assumed  to  be  in  a 
lightly  loaded  condition,  as  mentioned  earlier  under  the  Discussion  of 
Purposes  and  Methods.  A  specific  static  stability  boundary  was  calculated 
with  the  simple  energy  analysis;  the  same  values  of  tip  loss  and  structural 
root  cutout  as  those  for  the  extended  Normal  Mode  Transient  Analysis  cal¬ 
culations  were  used.  Ibis  is  the  boundary  that  appears  on  Figure  21. 

Almost  the  same  boundary  is  obtained  from  Figure  2  by  defining  the  stiffness 
coefficient  S*  in  terms  of  the  blade  length  from  structural  cutout  to 
the  outermost  aerodynamically  effective  radius.  When  the  strictly  uniform 
blade  without  tip  loss  or  structural  and  aerodynamic  root  cutout  is  con¬ 
sidered,  the  static  divergence  analysis  predicts  a  critical  center-of- 
gravity  offset  of  T%  chord.  This  offset  is  measured  aft  of  the  25/t  chord, 
where  the  advancing  blade  center  of  pressure  and  elastic  axis  are  located. 
The  boundary  shown  on  Figure  21  appears  at  a  center-of -gravity  offset  of 
10. 3£  chord.  Thus  a  reasonably  conservative  result  for  critical  center-of- 
gravity  location  can  be  obtained  by  ignoring  tip  loss  and  the  normal  amount 
of  structural  cutout. 
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The  information  on  Figures  22  and  23  is  for  the  nonarticulated  rotor.  It 
can  be  3een  that  the  torsional  stress  response  is  very  similar  to  that  for 
the  articulated  blade.  The  absence  of  the  blade  hinges  does,  of  course, 
cause  the  blade  bending  stresses  to  be  somewhat  higher.  It  should  be 
remembered  that  this  particular  result  is  for  a  uniform  blade  and  that 
design  refinements  can  alter  the  stress  distribution.  The  stress  results 
presented  here  are  intended  only  to  show  that  a  practical  operating  limit 
is  being  approached.  The  static  torsional  stability  boundaries  are  vir¬ 
tually  identical  for  articulated  and  nonarticulated  blades  at  p  ■  .6,  as 
shown  by  Figure  13. 

Figure  2k  shows  the  results  of  the  gyro  feedback  calculations.  It  can  be 
seen  that  a  fairly  large  increase  in  blade  torsion  response  occurs  because 
of  the  gyro  action  with  the  particular  control  gyro  parameters  chosen. 


DISCUSS IOWS  OF  TORSIONAL  PI V EROENCE  RESULTS 


SIGNIFICANCE  OF  TORSIONAL  DIVERGENCE  DESIGN  PARAMETERS 


The  scope  of  this  study  included  independent  variations  in  all  of  the 
important  parameters  affecting  the  torsionai  divergence  calculations  for 
helicopter  blades.  The  various  charts  included  in  this  volume  show  the 
relative  significance  of  the  various  design  parameters  as  well  as  quantita¬ 
tive  data.  The  practical  significance  of  changes  in  some  of  these  par¬ 
ameters  will,  however,  be  clarified  by  additional  discussion. 

The  basic  relationship  for  the  torsional  divergence  boundaries  is  between 
the  familiar  advance  ratio  fx  and  the  nondimens ional  stiffness  coefficient 
Sn  ,  which  is  basically  a  ratio  of  elastic  and  aerodynamic  torsional 
stiffness  effects.  It  should  be  noted  that  the  rotational  tip  speed  is 
used  in  the  formation  of  these  parameters.  Thus,  an  increase  in  Sr  can 
signify  a  larger  torsional  elastic  stiffness  or  a  decrease  in  the  quantity 
in  the  denominator  in  Eq.  (108).  This  may  occur  for  a  given  blade  because 
of  operation  with  lower  air  density  or  rotational  speeds.  As  would  obvi¬ 
ously  be  expected,  increasing  the  value  of  Sn  raises  the  critical  advance 
ratio  in  all  cases.  An  increase  in  S  „  resulting  from  a  decrease  in  ro¬ 
tational  speed  does,  however,  result  in  a  drop  in  critical  aircri.ft  forward 
speed  for  retreating  blade  divergence.  Therefore,  the  typical  appearance 
of  the  divergence  boundaries  of  a  representative  blade  on  a  tip  speed  versus 
forward  speed  plane  is  of  important  practical  interest.  Examples  of  this 
relationship  are  shown  in  Figures  25  through  28.  Tnese  charts  were  pre¬ 
pared  from  Figures  2,  3,  and  U,  with  Sn  assumed  to  be0.oUo  when  tip  speed 
is  660  ft/sec.  This  is  a  uniform  blade  comparable  to  the  CH-3C  blade  when 
operating  under  normal  conditions  with  a  lift-curve  slope  of  2  tt  . 

The  effect  of  changing  the  parameter  S*(RiI)2=  2GJo/padCo2R2  is  shown  on 
Figure  25.  Note  that  this  parameter  includes  blade  torsional  stiffness, 
chord,  radius,  air  density,  and  lift-curve  slope.  The  boundaries  shown  on 
the  chart  for  various  values  of  this  parameter  include  that  for  the  repre¬ 
sentative  blade  and  those  for  2/3  and  3/2  of  that  value.  Note  also  that 
any  number  of  these  boundaries  may  be  plotted  from  a  single  Sr  versus 
/i.  plot.  Thus,  this  choice  of  parameters  greatly  condensed  the  number  of 
charts  necessary  to  evaluate  the  effects  of  any  appreciable  variations  in 
2GJq/^>QqCq  R 

The  effect  of  blade  chordwise  center-of-gravity  position  is  shown  by  the 
various  corresponding  boundaries  in  Figure  26.  The  aft  center-of-gravity 
positions  give  rise  to  advancing  blade  boundary  segments,  which  appear  as 
lines  nearly  parallel  to  the  constant  advancing  tip  Mach  number  line.  The 
effect  of  changes  in  blade  lift-curve  slope  due  to  Mach  number  is  not  in¬ 
cluded  on  this  chart,  since  only  approximate  quantitative  information  is 
desired.  If  the  correction  for  Mach  number  is  desired,  the  actual  lift- 
curve  slope  for  a  given  Mach  number  at  the  75/J  radius  can  be  used  in  the 
definition  of  Sr  .  The  Mach  number  is  generally  low  on  the  retreating 
blade,  so  the  lift  curve  slope  on  the  corresponding  segment  of  the  boundary 
will  not  be  appreciably  affected  by  compressibility  effects. 
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The  moat  inboard  nondimensional  radius  for  which  significant  aerodynamic 
blade  torques  are  produced  is  called  the  aerodynamic  root  cutout,  X0A 
This  parameter  has  a  significant  effect  on  retreating  blade  torsional 
divergence,  since  th  inboard  blade  sections  are  subjected  to  the  highest 
velocity  reversed  flow.  It  can  be  seen  from  Figures  2,  3,  and  U  that  the 
beneficial  effect  of  root  cutout  is  more  pronounced  with  the  larger  values 
of  l/K,  as  would  air  >  be  expected,  since  a  smaller  control  stiffness  causes 
greater  blade  deflections,  particularly  close  to  the  root.  For  the  repre¬ 
sentative  blades.  Figure  27  shows  the  effects  of  changing  roc-  cutout  on 
the  flight  condition  boundaries. 

The  control  stiffness  ratio,  K:KfLB/J0G  ,  can  be  also  defined  as  blade 
deflection  divided  oy  the  control  system  deflection  when  the  blade  and 
control  system  are  subjected  to  a  torsional  couple  at  the  tip.  It  will  be 
noted  that  this  parameter  is  treated  in  inverse  form  on  the  charts  to  per¬ 
mit  plotting  very  large  values  of  K.  The  value  K  =  1,000  is  large  enough 
so  that  only  infinitely  small  effects  on  torsional  divergence  due  to  root 
flexibility  occur.  For  the  representative  blades,  Figure  28  shows  the 
effects  of  changing  control  stiffness  ratio. 

It  will  be  noted  from  Figures  25  through  28  that  the  practical  variations 
considered  for  root  cutout,  center-of-gravity  position,  and  control  stiff¬ 
ness  ratio  generally  have  only  moderate  effects  on  the  flight  condition 
boundaries  for  the  retreating  blade.  The  aft  center-of-gravity  positions 
do,  of  course,  make  an  advancing  blade  torsional  divergence  possible  at 
much  lower  speeds. 

The  ratio  of  blade  chord  to  rotor  radius  C0/R»  has  a  negligible  effect  on 
the  tvo-degree-of-freedom  torsional  divergence  analysis  for  variations  of 
C0/R  between  Q 0^  and  0. 06. 

The  independent  effect  of  linear  taper  in  chord  is  shown  in  Figure  5-  It 
should  be  remembered  that  these  variations  were  made  with  blade  chord  at 
the  75#  radius  remaining  constant,  and  with  the  torsional  stiffness  coef¬ 
ficient  defined  using  that  chord.  All  other  blade  properties  were  the  same 
as  for  the  corresponding  uniform  blade.  It  is  interesting  to  note  that 
with  the  above  definition  of  Sr  ,  the  addition  of  chord  taper  causes  a 
moderate  reduction  in  retreating  blade  divergence  advance  ratio.  However, 
if  Sr  is  efined  further  inboard,  the  boundaries  could  practically  coin¬ 
cide  or  ev  n  be  reversed  in  their  relative  positions.  The  practical  con¬ 
clusion  for  this  type  of  taper  is  better  suggested  by  Figure  5,  with  a 
moderate  reduction  in  divergence  speed  because  of  larger  inboard  chord  in 
the  region  of  higher  reversed  velocity. 

The  structural  cutout,  or  torsionally  rigid  inboard  section  has  the 
straightforward  effect  shown  in  Figure  6. 

The  independent  effect  of  tapered  tcrrionai  stiffness  on  retreating  blade 
divergence  is  shown  in  Figure  7-  It  can  be  seen  that  the  nondimensional 
boundary  is  virtually  unaffected,  provided  the  torsional  stiffness  coef¬ 
ficient  SR  is  defined  using  blade  stiffness  at  the  75#  radius  station. 
This  result  arises  because  sm-ller  stiffness  outboard  of  the  75#  radius 
counteracts  the  benefit  of  increased  stiffness  inboard. 
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The  independent  effects  of  nonuniform  inertia,  which  enter  into  the  tor¬ 
sional  divergence  static  stability  calculations  because  of  variations  in 
the  assumed  torsional  mode  shape,  are  shown  in  Figure  8  and  9*  Figure  8 
was  prepared  for  two  values  of  aerodynamic  cutout  by  assuming  an  inboard 
section  with  negligible  torsional  Inertia.  It  can  be  seen  that  this  vari¬ 
ation  has  only  a  small  effect.  As  would  be  expected,  the  generally  smaller 
deflections  inboard  for  the  mode  shape  of  the  blade  with  inertir  cutout 
cause  a  higher  critical  advance  ratio  to  be  predicted. 

The  effect  of  a  linear  variation  in  inertia  is  shown  in  Figure  9 .  This 
figure  was  prepared  by  assuming  torsional  inertia  at  the  75%  radius  equal 
to  that  of  the  uniform  blade.  It  can  be  seen  that  the  generally  larger 
inboard  deflections  for  the  mode  shape  of  the  blade  with  inertial  taper 
cause  a  small  decrease  in  predicted  divergence  advance  ratio. 

Spanwise  variations  in  center-of-gravity  location,  with  the  location  at  the' 
75%  radius  remaining  constant  at  the  10£  chord  forward  position,  were  con¬ 
sidered.  Inboard  cutouts  of  center-of-gravity  offset  up  to  the  b5%  radius 
and  linear  tapers  with  the  ratio  of  root  offset  to  tip  offset  as  high  as 
three  produced  virtually  no  effect  on  the  retreating  blade  cases  considered. 

The  effect  on  advancing  blade  divergence  of  a  lift-curve  slope  distribution 
due  to  the  familiar  Glauert  Mach  number  correction  is  shown  in  Figure  10. 
Vote  that  the  torsional  stiffness  coefficient  is  defined  by  using  the  local 
lift-curve  slope  at  the  75%  radius.  The  value  Sc  ■  .0286  refers  to  a 
blade  similar  to  a  conventional  blade  with  twice  the  torsional  stiffness 
operating  under  normal  conditions.  It  can  be  seen  that  with  torsional 
stiffness  coefficient  defined  in  the  manner  described  earlier,  very  little 
effect  due  to  Mach  number  lift-curve  slope  variation  can  be  noted  on  the 
nondlmenslonal  charts.  Even  less  effect  is  noted  on  the  retreating  blade, 
shown  in  Figure  11,  as  would  be  expected.  Even  though  the  radial  variation 
in  Mach  number  effects  has  little  effect  on  the  nondlmenslonal  charts,  it 
should  be  remembered  that  an  effect  will  be  felt  on  the  actual  flight  con¬ 
dition  boundary,  since  a  larger  average  lift-curve  slope  is  reflected  in  a 
smaller  stiffness  coefficient. 

The  effect  of  bending  mode  shape  for  a  nonart iculated  blade,  as  compared 
with  the  rjid  flapping  mode  used  for  the  articulated  blade,  is  shown  in 
Figure  12  and  13.  It  should  be  noted  that  the  uniform  blade  natural  bend¬ 
ing  mode  shape  and  natural  frequency  to  rotational  speed  ratio  depend  only 
on  the  parameter  EI^/mH^R  ,  so  that  a  definite  mode  shape  and 

frequency  ratio  is  defined  by  the  product  of  S  «  ,  the  ratio  of  bending 

to  torsion  stiffness  El, /GJ0  *  and  a  parameter  p  o0c0z/m  .  The 
effect  of  articulation  was  determined  for  tw6  values  of  EIf/GJ0  with  the 
standard  value  of  p o0c02/m  .  The  same  variations  in  critical  advance 
ratio  would  occur  for  a  change  in  the  product  of  these  two  parameters 
EIrpOoCo*/GJow  by  the  same  factors.  The  value  EIp/GJ0  of 0. 598  Is 

similar  to  conventional  practice,  while  p o0t£tm  is  also  appropriate  to 
a  lift-curve  slope  of  2  n  and  standard  sea  lev*l  conditions.  It  is  readily 
apparent  that  the  values  considered  produce  very  little  effect  on  torsional 
divergence  for  the  retreating  blade  and  a  moderate  effect  for  the  advancing 
blade.  The  changes  in  the  relative  positions  of  the  various  boundaries  are 


probably  caused  by  interacting  favorable  effects  of  bending  stiffness  and 
unfavorable  effects  of  the  various  bending  mode  shapes. 

The  pronounced  lowering  of  retreating  blade  torsional  divergence  advance 
ratio  predicted  by  the  simple  energy  method  for  pitch-flap  coupling  is 
shown  in  Figure  lh.  It  should  be  noted  that  the  studies  with  the  extended 
Normal  Mode  Transient  Analysis  show  this  predicted  effect  to  be  very  on- 
servative.  The  effects  of  pitch-flap  coupling  for  the  initially  unloaded 
rotor  can  be  seen  by  comparing  Figures  17  and  20.  The  peak-to-peak  tor¬ 
sional  stress  is  actually  reduced  by  pitch-flap  coupling  for  a  given  advance 
ratio  while  the  peak-to-peak  flatwise  stress  increases.  However,  the  practical 
retreating  blade  advance  ratio  stress  boundary  is  very  little  affected  by 
the  addition  of  Sj  *  U50 ,  according  to  these  Normal  Mode  Transient 
Analysis  results. 

The  provision  of  the  control  gyro  configuration  used  in  this  study  was 
found  to  increase  blade  torsional  response,  and,  therefore,  to  lower 
slightly  the  permissible  advance  ratio  for  a  given  stress  by  a  decrement  of 
approximately  .10.  This  is  shown  in  Figures  22  and  2h.  It  should  be  remem¬ 
bered  that  the  effect  of  the  gyro  can  be  changed  by  altering  the  various 
parameters  and  that  no  general  implications  can  be  drawn  without  further 
analyses . 

It  has  been  demonstrated,  as  discust-ed  earlier  in  this  section,  that  the 
most  significant  improvements  in  torsional  divergence  speed  may  be  brought 
about  by  changing  the  quantities  involved  in  the  torsional  stiffness  coef¬ 
ficient  Ss  .  Thus,  an  increase  in  the  quantity  GJb//°°oco,R  ,  which  can 
be  brought  about  by  an  increase  in  blade  uorsional  stiffness  or  a  decrease 
in  chord,  rotor  radius,  lift-curve  slope,  01  air  density,  will  raise  the 
speed.  On  the  other  hand,  decreasing  rotor  rotation'-l  speed  will  lower  the 
aircraft  divergence  speed  for  retreating  blade  divergence,  but  will  increase 
it  for  advancing  blade  divergence.  If  the  blade  has  an  excessively  large 
aft  chordwise  center-of-gravity  position,  very  large  improvements  in  tor¬ 
sional  divergence  speed  are  possible  by  returning  to  the  normal  balanced 
blade  configuration.  Significant  improvements  may  also  b«*  realised  by  pro¬ 
viding  stiff  inboard  blade  sections,  as  illustrated  by  Figure  6.  This 
should  not  be  undone  by  decreasing  the  stiffness  of  outboard  blade  sections. 

In  Figure  f ,  it  can  be  seen  that  the  divergence  boundaries  for  a  tapered 
blade  and  a  uniform  blade  are  almost  identical.  The  beneficial  high  inboard 
stlffneaa  Indicated  by  Figure  6  must  be  derived  from  an  Inboard  atlffnass  that 
is  higher  in  an  absolute  sense  rather  than  higher  relative  to  the  tip.  1c 
should  he  remembered  that  Figure  ^  was  prepared  with  the  aame  atlffness  at  the 
75$  radius  for  tapered  and  untapered  blades,  with  all  other  parameters,  includ¬ 
ing  chord,  remaining  unchanged.  Other  parameters  investigated  ware  found  to 
have  moderate  Independent  effects,  but  it  can  be  expected  that  significant  im¬ 
provements  will  be  achieved  through  the  combined  favorable  effects  of  changes 
in  a  nuafcer  of  them.  The  significance  of  the  various  parameters  investigated 
has  been  discussed  above,  with  respect  to  the  torsional  divergence  calcu¬ 
lated  results.  The  practical  significance  of  the  torsional  divergence 
concept  itself  should  be  reviewed,  as  it  applies  to  a  helicopter  rotor 
blade.  The  most  superficial  inspection  of  the  Normal  Mode  Transient  Anal¬ 
ysis  results,  such  as  Figures  15  and  16,  sh^ws  that  the  calculated  stability 


boundary  for  torsional  divergence  corresponds  generally  to  the  onset  of 
torsional  oscillations  of  comparatively  high  frequency.  Thus,  the  basic 
equations  of  static  torsional  divergence,  which  neglect  all  inertial 
effects,  cannot  be  viewed  as  a  mathematical  model  of  a  physical  phenomenon 
in  the  usual  sense.  The  torsional  divergence  calculated  results  should 
rather  be  viewed  as  a  so-called  "rule  of  thumb".  The  virtue  of  these  re¬ 
sults  lies  in  the  relatively  simple  manner  in  which  they  can  be  obtained 
and  the  fact  that  the  results  do  indeed  correspond  with  those  of  a  much 
more  elaborate  method.  The  practical  boundary  for  retreating  blades,  as 
limited  by  peak-to-peak  torsional  stress,  is  accurately  predicted  for  all 
parameter  variations  considered  with  the  Normal  Mode  Transient  Analysis, 
except  for  the  effect  of  pitch- flap  coupling.  The  aft  center-of-gravity 
advancing  blade  stress  "boundary"  is  not  as  well  defined,  probably  because 
of  complicated  compressibility  effects.  The  torsional  static  stability 
boundary  provides  a  conservative  estimate  of  the  rather  gentle  onset  of 
high  torsional  stresses  for  the  unloaded  blade. 

Therefore,  in  view  of  the  above,  the  inclusion  cf  a  large  number  of  elab¬ 
orate  features  in  a  torsional  divergence  calculation  is  apparently  not 
warranted.  It  should  also  be  pointed  out  that  predicted  trends  for  par¬ 
ameters  not  investigated  previously  may  not  necessarily  be  accurate  or 
correct.  Torsional  divergence  results  for  parameter  variations  not  yet 
considered  and  for  radically  different  designs  should  therefore  be  checked 
with  more  elaborate  methods,  as  was  done  under  this  contract  for  some  par¬ 
ameters.  The  effect  of  inertia  distributions  on  torsional  divergence  should 
be  checked  with  particular  thoroughness. 

The  ultimate  objective  of  torsional  divergence  calculations  should  always 
be  a  guide  for  more  elaborate  methods .  When  a  relatively  large  number  of 
torsional  divergence  calculations  have  defined  flight  conditions  and  con¬ 
figurations  critical  for  blade  torsional  stiffness,  a  much  smaller  number 
of  dynamic  solutions  can  be  used  to  check  these  critical  regions. 

PERMISSIBLE  ANALYTICAL  ASSUMPTIONS 


In  the  following,  the  torsional  divergence  calculations  and  results  will  be 
reviewed  with  the  objective  of  outlining  useful  analytical  assumptions  for 
such  calculations. 

The  results  shown  in  Figures  2  through  U  establish  that  the  basic  advance 
ratio  versus  torsional  stiffness  coefficient  relationship  is  affected  by 
practical  variations  in  the  aerodynamic  root  cutout,  control  stiffness 
ratio,  and  blade  chordwise  center-of-gravity  position.  The  effects  of  these 
parameters  are  interrelated  in  such  a  way  that  no  generally  useful  statement 
concerning  negligible  values,  which  will  be  suitable  for  all  cases  can  be 
given.  The  effect  of  root  cutout,  for  example,  is  more  pronounced  when  the 
control  stiffness  ratio  is  low. 

Figures  $  and  6  show  that  spanvise  variations  in  chord  and  structural  cut¬ 
out  have  moderately  Important  effects. 
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Figure  7,  on  the  other  hand,  shows  that  linear  variations  in  torsional 
stiffness  msy  be  neglected  in  the  simple  torsional  divergence  calculation, 
if  the  value  of  stiffs 3ss  at  the  75%  radius  is  used  to  represent  an  equiv¬ 
alent  blade  with  unifi  nn  stiffness. 

Figures  8  and  9  reflect  variations  in  torsional  mode  shape  consistent  with 
quite  extensive  variations  in  blade  mass  distributions.  These  can  be  seen 
to  produce  rather  small  changes  in  calculated  torsional  divergence  results, 
and  the  use  of  a  uniform  blade  torsional  mode  shape  will  be  adequate  for 
most  cases.  This  would  not  necessarily  be  true  for  even  more  extensive 
variations  in  modal  properties,  such  as  those  appropriate  to  a  weightless 
blade  with  tip  mass.  It  should  be  pointed  out  again  that  the  effect  of 
mass  on  torsional  divergence  is  not  logically  consistent  with  the  basic 
assumption  of  no  inertial  effects.  Earlier  work,  however,  has  indicated 
that  the  correct  trends  appear  to  be  produced  by  using  the  natural  mode 
shape  corresponding  to  the  mass  distribution  in  the  static  torsional  diver¬ 
gence  calculations.  While  a  basic  inconsistency  is  present  from  a  strictly 
logical  point  of  view,  the  use  of  torsional  divergence  calculations  for 
nxle-of- thumb  purposes  makes  this  practice  reasonable,  since  the  effect  of 
mass  distribution  on  torsional  stress  appears  to  be  predicted  with  the 
desired  accuracy.  This  characteristics  of  the  torsional  divergence  calcu¬ 
lations  should  be  checked  with  more  elaborate  analyses,  as  mentioned  pre¬ 
viously. 

Figures  10  and  11  show  that  lift-curve  slope  corrections  may  be  confined  to 
the  lift-curve  slope  at  the  75%  radius  in  the  definition  of  torsional  stiff¬ 
ness  coefficient  Sr  .  This  will  be  true  if  the  spanwise  variation  is  no 
more  severe  than  that  due  to  the  Prandtl-Glauert  correction  for  an  advancing 
blade  75%  radius  Mach  number  of  approximately  0.70. 

Figures  12  and  13  show  that  nonarticulated  blades  of  conventional  design 
will  have  calculated  torsional  divergence  speeds  close  to  those  for  artic¬ 
ulated  blades.  This  is  confirmed  by  the  Normal  Mode  Transient  Analyses 
results  shown  in  Figures  17,  21,  22,  and  23. 

Figure  1U  shows  a  large  drop  in  torsional  divergence  advance  ratio  due  to 
pitch- flap  coupling,  which  was  not  reflected  in  the  results  of  the  Normal 
Mode  Transient  Analysis  shown  in  Figures  17  and  20.  The  effects  of  85 
should  therefore  be  investigated  with  a  more  rigorous  analysis. 

In  addition  to  the  above  considerations,  the  independent  effects  of  varia¬ 
tions  in  mass  ratio  Mr  were  found  to  be  negligible.  The  effects  of 
linear  spanwise  variation  in  chordwise  center-of-gravity  position  were  also 
small,  if  the  position  at  the  75%  radius  was  used  to  compare  with  the  uni¬ 
form  blade  results. 

TORSIONAL  DIVERGENCE  PRELIMINARY  DESIGN  PROCEDURES 


The  following  paragraphs  include  recommendations  as  to  the  use  of  the  tor¬ 
sional  divergence  charts  in  the  preliminary  design  stages  of  new  rotor 
systems . 
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It  should  be  recognized  that  a  complete  outline  of  preliminary  design  pro¬ 
cedure  is  beyond  the  scope  of  this  study.  Examples  of  hov  information 
useful  in  preliminary  design  may  be  obtained  from  the  charts  are  given. 

The  designer  should  determine  the  specific  studies  to  be  made.  These 
studies  must  determine  advantageous  combinations  including  not  only  required 
stiffness  but  also  items  such  as  weight,  cost,  and  aerodynamic  performance. 

The  torsional  divergence  charts  can  be  used  to  define  allowable  relation¬ 
ships  between  torsional  stiffness,  blade  geometry  and  operating  conditions. 
It  can  be  expected  that  large  torsional  stresses  will  arise  if  these  re¬ 
lationships  are  exceeded  significantly. 

It  should  be  remembered  that  high  torsional  stresses  may  arise  for  other 
reasons,  even  though  the  blade  is  satisfactory  from  the  standpoint  of  tor¬ 
sional  divergence.  The  presence  of  classical  or  stall  flutter,  vhich  may 
occur  well  below  the  torsional  divergence  boundary,  will  also  result  in 
high  torsional  stresses  .  Excessive  flapping  motion  may  also  result  in 
large  torsional  stresses,  depending  on  the  modal  coupling  present.  Tor¬ 
sional  resonance  with  rotor  orders  must  also  be  avoided. 

It  must  be  assumed  that  certain  parameters  have  been  defined  at  the  start 
of  a  preliminary  design  effort  by  the  basic  mission  requirements  and  the 
resulting  aerodynamic  preliminary  rotor  design.  In  particular,  it  is 
assumed  that  rotor  diameter,  blade  chord,  aircraft  forward  velocity,  rotor 
rotational  speed,  and  air  density  are  known  or  chosen.  In  preliminary 
design,  it  is  recommended  that  Figures  2  through  k  be  used.  If  the  normal 
amounts  of  tip  loss  and  structural  root  cutout  are  neglected,  an  appropriate 
amount  of  conservatism  will  be  introduced. 

With  the  quantities  .Just  mentioned,  the  advance  ratio  is  calculated,  and 
Figures  2,  3»  or  U  entered  to  find  a  required  Sr  for  the  desired  values 
of  aerodynamic  zoot  cutout  XM  .  Then  the  definition  of  Sr  given  in 
Eq.  (108)  is  used  to  find  the  torsional  stiffness  GJ0  required.  The 
lift-curve  slope  o<>  is  appropriate  for  the  conditions  at  the  75f  blade 
radius  on  the  advancing  or  retreating  blade. 

The  above  procedure  can  be  carried  out  quickly  for  systematic  variations 
in  control  stiffness  ratio  and  center-of-gravity  location  by  using  the 
trade-off  charts  in  Figure  3  and  l*.  The  blade  weight  for  each  combination 
of  torsional  stiffness,  control  stiffness  ratio,  and  center-of-gravity 
position  can  then  be  determined  to  find  the  minimum  weight  configuration 
for  the  given  flight  condition. 

Another  possible  procedure  might  begin  with  blade  stiffness,  chord,  approx¬ 
imate  lift-curve  slope,  and  radius  defined,  so  that  variations  in  the 
quantity  Sr  (  Rfl  )*  are  determined  by  air  density  and  rotational  tip 
speed.  Then  flight  condition  boundaries  similar  to  those  shown  in  Figures 
2$  through  28  can  be  determined  for  various  values  of  altitude,  root  cut¬ 
out,  or  control  stiffness  ratios.  This  is  done  by  choosing  a  list  of  ad¬ 
vance  itlos,  determining  the  corresponding  list  of  critical  Sr  values 
from  Figures  2,  3,  or  l,  and  then  solving  for  rotational  tip  speed  HR 
from  the  value  of  Sr(RA)*  ,  consistent  with  the  assumed  altitude.  Mote 


that  the  altitude  is  reflected  in  the  value  of  air  density.  The  exact  Mach 
number  and  lift-curve  slope  will  not  be  known  until  the  boundary  is  deter¬ 
mined  if  this  method  is  used;  therefore,  it  may  be  desired  to  repeat  the 
process  with  the  lift-curve  slope  appropriate  to  the  approximate  Mach 
number  for  each  of  the  advance  ratio  and  rotational  tip  speed  solutions. 

An  alternate  procedure  to  obtain  flight  condition  boundaries,  such  as  those 
described  ir  the  preceding  discussions,  is  based  on  the  assumption  of  vari¬ 
ous  advancing  blade  tip  Mach  numbers.  This  procedure  permits  the  use  of 
the  correct  lift-curve  slope  without  need  for  trial  and  error,  but  it  is 
slightly  more  elaborate.  This  method  is  especially  useful  for  the  high 
advance  ratio  and  high  Mach  number  conditions. 

The  advancing  blade  tip  and  75%  radius  Mach  numbers  on  the  advancing  and 
retreating  blades  sure  given  respectively  by 


M 


i.*o 


flR(l+M) 

A* 


(138) 


M.rs.to 

M.ts.ito 


flR(.75+M)  ..  (75*m> 

A«  •  m'.“TTt£t 

GRIM"75)  ..  I^-75| 

A*  ‘  '•#0  I4-/4 


(139) 

(lUo) 


When  Eq.  (138)  is  used  in  the  definition  of  SR, 

S  -  2GJ«  .  ZSJaU+M)* 

"  "  ,a0C0iRi(flR}1  '  ,OoCjR'(AX»o> 


(1«*1) 


When  the  advancing  blade  tip  Mach  number  Ml#0  and  the  value  of  advance 
ratio  pi  are  defined,  Eq.  (lbl)  can  be  used'  to  calculate  a  value  of  Sa  . 
The  lift-curve  slope  at  the  75%  radius  is  used  in  this  calculation.  The 
lift-curve  slope  is  corrected  for  the  Mach  number  at  the  75%  radius  which 
is  calculated  by  using  Eq.  (139)  or  Eq.  (lbo)  for  the  advancing  or  retreat¬ 
ing  blade  respectively.  Thus,  two  curves  of  pi  versus  S»  are  obtained, 
one  for  the  retreating  blade  and  one  for  the  advancing  blade  if  there  is 
an  aft  center-of-gravity  location.  These  are  plotted  on  the  appropriate 
chart  in  Figure  2,  and  the  intersection  with  the  divergence  boundaries  is 
noted.  The  same  pi  and  intersection  is  then  plotted  on  a  rota¬ 

tional  tip  speed  flR  versus  forward  speed  iV)  chart  to  give  a  point  on 
the  flight  condition  boundary.  This  is  repeated  for  other  values  of  M(t0 
until  a  torsional  divergence  boundary  is  defined.  Other  boundaries  can' be 
determined  as  desired  by  using  a  new  value  of  the  quantity  2GJ0/pCe*  R* 
in  Eq.  ( lhl ) .  Obviously,  changes  in  this  group  of  parameters  can  come 
about  because  of  variations  in  torsional  stiffneso,  air  density,  and  blade 
area.  The  effects  of  aerodynamic  root  cutout,  control  stiffness  ratio,  and 
center  of  gravity  are  considered  by  using  the  proper  curves  in  Figure  2. 


Approximate  corrections  for  chcrd  taper,  structural  cutout.  Inertia  varia¬ 
tions  and  nonarticulation  may  be  made  if  desired.  This  can  be  done  by 
adding  the  appropriate  increment  to  Se  ,  as  indicated  at  each  value  of 
fi  in  Figures  5  through  10  and  12  and  13. 

A  few  procedures  have  been  outlined  above  for  the  rapid  definition  of  sat¬ 
isfactory  relationships  between  blade  torsional  stiffness  and  rotor  per¬ 
formance  objectives  using  the  design  charts  presented.  These  should  not 
be  considered  as  restrictions,  however,  since  other  procedures  will  probably 
occur  to  various  users.  Detailed  explanations  of  how  the  charts  were  ob¬ 
tained  have  been  given  to  allow  individual  users  to  develop  their  own  pro¬ 
cedures.  The  major  restriction  on  the  development  of  other  procedures  is 
that  they  be  simple,  rapid,  and  consistent  with  the  approximate  nature  of 
the  torsional  divergence  concept  for  helicopter  blades.  In  the  later  stages 
of  design  work,  the  more  critical  flight  conditions  must  be  reevaluated 
with  more  elaborate  methods,  including  those  for  the  prediction  of  forced 
blade  motions  and  stresses,  oxassical  flutter,  and  stall  flutter. 
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Figure  1.  Comparison  of  Differential  Equation  Solution  and  Assumed 
Divergence  Mode  . 
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Figure  2.  Advance  Ratio  Versus  Torsional  Stiffness  Coefficient  fcr 
Various  Control  Stiffness  Ratios  and  Center-of-Gravity 
Positions;  XM  =  0  ,  o0  -  Zir  . 
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FIGURE  2.  CONTINUED. 
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FIGURE  2.  CONTINUED. 
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FIGURE  2.  COUTIKUED. 
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FIGURE  2.  COIITINUED. 
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FIGURE  2.  CONTINUED. 
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Figure  5.  Effect  of  Linear  Variation  in  Chord  on  Retreating  Blade 
Toruional  Divergence  ;  l/K  *  001  ,  YcG^C  =  0  »  xo*  :  25  , 

Xp»  :  0  ,  00  :  2  7T  . 
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Figure  «  . 


Effect  of  Structural  Cutout  on  Hetreatinp  Blade  Torsional 
Divergence;  l/K  :  .001  ,  Ycc/C  =0,  Xe*  =  0 ,  o0  =  2  tt  . 
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Figure  7.  Effect  of  Linear  Variation  in  Stiffness  on  Retreating 
Blade  Torsional  Divergence}  l/K  =  .001  ,  Ycg/C  5  0  , 

*oa  :  0  ,  X0,  =  0  ,  o0  =  2tt  . 
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Figure  8.  Effect  of  Inertial  Cutout  on  Retreating  Blade  Torsional 
Divergence,  l/K  -  001  ,  Ycg/C  =  0  ,  X05  =  0  ,  o0  -  Zv  . 
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Figure  9.  Effect  of  Linear  Variation  in  Inertia  on  Retreating 
Blade  Torsional  Divergence;  l/K  =  .001  ,  X0A  =  0  , 

Xqj  -  0  ,  Oo  •  2  7T 
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Figure  10.  Effect  of  Franitl-Glauert  Mach  Number  Correction  on 
Advancing  Elade  Torsional  Divergence;  l/K  =  .OOI  , 
Yqg/C  =  —  . 1 0  ,  XqazO,  Xos  =  0 
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Figure  11.  Effect  of  Prandtl-Glauert  Mach  Number  Correction  on 
Retreating  Blade  Torsional  Divergence;  I / K  =  .001  , 
Yro/C  -  0  ,  X0A  =  0  ,  Xos  =  0 
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Figure  12.  Retreating  Blade  Torsional  Divergence  on  Articulated  and 
Nonarticulated  Blades;  l/K  =  .001  ,  Ycg/C  =  -005  ,  X0A  =  0  , 
Xq$  “0  ,  p  o0C0*/  m  =  .168  . 
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Figure  13.  Advancing  Blade  Torsional  Divergence  on  Articulated  and 
Nonarticulated  Blades  j  l/K  =.001  ,  Ycg/C=  -.005  ,  XM  =0  , 
Xos  =  0  ,  p  OoCo*/m  =  .168  . 
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Figure  15*  Blade  First  Torsion  Modal  Response  Versus  Azimuth  Angle  for 
the  Basic  Articulated  Rotor;  ft  =  1.6  ,  d,7in  :  0°  ,  a  =0, 
Ycg / C  =  0  ,  M|,jo  :  85  .  ® 
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Figure  16.  Blade  First  Torsion  Modal  Response  Versus  Azimuth  Angle  for 
the  Articulated  Rotor  with  Aft  Center  of  Gravity;  fj.  -  .6  , 
0.75R  =  2°  ,  °s  =  0,  YCG/C  =  -.15  ,  Ml(,o  =  .85 
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Figure  17.  Vibratory  Stress  Versus  Advance  Ratio  for  the  Basic 
Articulated  Rotor;  l/K  =  0  ,  Ycg/C  =0  ,  X0A  =  .12  , 

X0s  "  .060  ,  M|,9o  —  .85  ,  S3  ~  0  ,  Ogil  Sr  -  96 
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Firure  18.  Vibratory  .'tress  Versus  Advance  natio  for  ‘he  Basic 
Articulated  Rctcr  with  y^edueei  Air  Density;  l/K  =  0  , 

Yc6/C  =0  ,  X0A  =  .12  ,  x0s=  060  ,  M|,*0  :  85  , 
p  =  001585  SLUGS /FT 3  ,  S3  =  0°  ,  o0fl*S*  =  |44. 
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Figure  19-  Vibratory  Stress  Versus  Advance  Ratio  for  the  Basic 
Articulated  Rotor  With  Control  System  Flexibility; 
l/K  =  2,  YCG/C  =  0  ,  X0A  =  12  ,  Xos  =  060  ,  M,t90  =  .85  , 

Sj  -  0°  ,  Coil  Sr  -  96  . 
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Vibratory  .'tress  Versus  Advance  Ratio  for  the  Basic 
Articulated  Rotor  With  Pitch-Flap  Courling;  l/K  =  0, 

Ycg/C  =0,  Xoa=  .12  ,  Xos  =  .060  ,  Mi>m  =  .  85  , 

Ss  =  45°  ,  a0fl‘SR  =  96  . 
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Figure  1  - . 


Vibratory  .'tress  Versus  Center— of- Iravi ty  OiTset  Fcr  the 

Art  i  culated  Rotor ;  l/K  =  0  ,  Xo*  =  12  ,  Xos  =  .060  ,  z  .  6  , 
Mi,»o  =  .85  ,  Sj  -  0  ,  o0il*S»  =  96  ,  079(,  =  2°  - 
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Figure  22.  Vibratory  ftress  Versus  Advance  Ratio  for  the  Basic 
I.'or.arti cuiatei  Rotor;  l/K  =  0,  Ycq/C  :  0  ,  X0A  =  12  , 

Xqs  ■  .060  |  ^i,9o  ”  85  f  o o  ii  Sr  =  96  , 
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Figure  23.  Vibratory  Stress  Versus  Center-of-Gravity  Offset  for  the 
Nonarticulated  Rotor,  I / K  =  0  ,  X 0A  =  .12  ,  Xos  =  060  , 
fi  =  6  ,  Ml<90  =  .85  ,  a0iIzSR  :  96  ,  075R  :  2°,  u  $  *  0,  Ais  " |C  , 
Bis  =~2V" 
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Figure  PC.  Effect-  of  Chordvise  Center  of  Gravity  Location  on 
Typical  Torsional  Divergence  Flight  Condition 
Boundaries ;  l/K  =  001  ,  X©a  *  25  ,  X©,  =  0  ,  o0  =  2tt  , 
S*(RQ)*  =  2GJ/^0oCo,R*  =  16,750  FT /SEC*  . 
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Figure  27.  Effects  of  Aeroivnanic  Root  Cutcut  on  Typical  Torsional 
Divergence  Flight  Condition  Boundaries;  l/K  =  .001  ,  Ycg/C  = 
Xos  =  0  ,  o0=  2tt  ,  S„{Rft)*  =  2GJ //>o0C02R*  =  16,750 FT/ SEC* 
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Figure  28.  Effects  of  Control  Stiffness  on  Typica^  Torsional  Divergence 
Flight  Condition  Boundaries;  YCg/C  =  0  ,  XM  =  .25  ,  X©*  =  O  , 

o0  =  2tt  ,  S„(Rfl)*  =  2GJ  / pOoZ*  R*  =  16 ,750  FT/SEC*  . 


TABLE  I.  BASIC  PROPERTIES  OF  ROTOR  SYSTEM  USED  IN  TORSIONAL  DIVERGENCE 
STUDIES  WITH  THE  NORMAL  MODE  TRANSIENT  ANALYSIS 


It«n 

Symbol 

Value 

Unit 

Radius 

R 

31.0 

ft 

Number  of  Blades 

b 

5 

- 

Chord 

c 

1.52 

ft 

Flatwise  Stiffness 

EIP 

1. 56x10* 

lb-ft* 

Edgewise  Stiffness 

1.7loloP 

lb-ft* 

Torsional  Stiffness 

GJ 

2.1*3x10* 

lb-ft* 

Hass  per  Unit  Length 

■ 

.202 

slugs/ft 

Flap-Lag  Hinge  Offset 

s’ 

.0339 

- 

Tip  Loes  Factor 

B 

.97 

- 

Aerodynamic  Root  Cutout 

*0. 

.12 

Structural  Root  Cutout 

lot 

.060 

TAELS  II.  ROTOR  CONFIGURATIONS  FOR  TORSIONAL  DIVERGENCE  STUDY  WITH 
THE  NORMAL  MODE  TRANSIENT  ANALYSIS 


Configuration 

lfuribar 

Nonartleulated 

1 

Articulated 

.001 

0.0 

0.0 

- 

2 

Artleulatad 

.20 

0.0 

0.0 

• 

3 

Articulated 

.001 

1.0 

0.0 

- 

h 

Articulated 

.001 

0.0 

-.05 

m 

$ 

Artleulatad 

.001 

0.0 

-.10 

m 

6 

Artleulatad 

.001 

0.0 

-.15 

- 

7 

Nonartlculatad 

.001 

- 

0.0 

fixed 

P. 

Nonartleulated 

.001 

- 

\Tv 

c 

• 

1 

Fixed 

9 

Nonartlculatad 

.001 

- 

-.10 

Fixed 

10 

Nonartleulated 

.001 

- 

-.15 

Fixed 

11 

Nonartleulated 

.001 

- 

0.0 

Acting 
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TABLE  III. 

FLIGHT  CONDITIONS  FOR  TORSIONAL  DIVERGENCE  STUDY 
WITH  NORMAL  MODE  TRANSIENT  ANALYSIS 

Condition 

Kinbor 

Adraneo  Ratio, P 

Air  D«nilty,  p 
Slugs/ft* 

1 

0.6 

.008378 

2 

uo 

•008376 

3 

1.2 

.008376 

li 

1.1 

.008376 

5 

1*6 

.002378 

6 

l.h 

.0015P5 

7 

1.6 

.001585 

8 

1.8 

.001585 
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CONCLUSIONS 


1.  The  static  torsional  divergence  boundary  usually  approximates  a 
practical  operating  limit  defined  by  a  rapid  rise  in  peak-to-peak 
torsional  stress.  This  approximation  is  valid  for  a  useful  number 
of  blade  parameters. 

2.  Critical  advance  ratio  for  static  torsional  divergence  is  affected 
nest  strongly  by  blade  elastic  torsional  stiffness,  air  density, 
lift-curve  slope,  rotational  tip  speed,  chord,  radius,  chordwise 
center-of-gravity  position,  aerodynamic  root  cutout,  control  stiff¬ 
ness,  and  structural  root  cutout. 

o.  The  torsional  divergence  characteristics  for  articulated  and  non- 
articulated  blades  of  conventional  stiffness  are  practically 
identical. 

**.  The  basic  parameters  for  the  study  of  static  torsional  divergence 
are  the  advance  ratio  ft  and  the  torsional  stiffness  coefficient 

S.  . 

5.  The  use  of  the  first  natural  vibration  mode  as  a  divergence  mode 

is  a  good  simplifying  assumption  for  static  divergence  calculations. 
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RECOMMENDATIONS 


The  present  static  stability  analysis  is  adequate  for  the  rapid  generation 
of  approximate  blade  operating  boundaries  for  a  broad  range  of  parameters 
and  blade  characteristics.  In  order  to  predict  the  effects  of  pitch-flap 
coupling  and  inertia  variations  with  greater  confidence,  a  tvo-degree-of- 
freedom  torsion-flapping  dynamic  analysis  should  be  considered. 
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^A  slaple  energy  balance  Method  for  datcnslnlng  static  stability  was  used  to  evsluete 
rotor  blade  torsional  divergence  boundaries,  and  a  set  of  blade  design  charts  was 
generated.  To  coopers  the  results  of  this  sleple  analysis  with  a  wore  refined 
method,  and  to  Investigate  the  effects  of  parameters  not  included  In  the  simple 
method,  other  calculations  were  made  with  the  extended  Normal  Mode  Transient  Analysis 

(hie  type  of  practical  rotor  blade  operating  boundary  Is  defined  by  peak-to-peak 
torsional  stress.  Ihe  static  torsional  divergence  stability  boundary  usually 
approximates  this  practical  operating  limit,  as  vsrifled  by  the  Normal  Mode  Transient 
Analysis  calculations.  This  approximation  is  valid  for  an  Important  number  of  blade 
paraawtera. 

Recommended  preliminary  design  procedures,  Incorporating  the  charts  generated  for 
this  volume  are  given  for  the  rapid  definition  of  blades  with  adequate  torsional 
stlffneas. 


